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The single particle stability in a circular accelerator is of concern especially for operational regimes 
involving beam storage of hours. In the proximity to a resonance this stability domain shrinks, and 
the phase space fragments into a jungle of exotic objects like for instance “fix-lines”. The concept 
of fix-points is easily understandable in a 2D phase space. It becomes quite challenging when the 
effect of resonances is considered in the 4D phase space, which leads then to the concept of fix-lines. 
In this paper we investigate the fix-lines in the proximity of a coupled third order resonance and 
find the relation of these objects with the stability of motion. 

PACS numbers: 41.75.-i, 29.27.Bd 


I. INTRODUCTION 

The stability of particle motion in circular accelera¬ 
tors has always been a topic of heated discussion in the 
accelerator community. In fact, the problem of localized 
magnet non-linearities distributed around the accelerator 
has invoked its own held of nonlinear dynamics studies. 
The most complex consequences of these non-linearities 
are detuning with amplitude, excitation of resonances, 
chaotic motion and eventually particle loss [HE]. 

The theory of resonances was developed by Hage- 
dorn/Schoch |5H5| in the 1950’s and modihed for CERN 
Accelerator School classes by Guignard in the 1970’s [SI 
[7] . This theory treats the condition under which the mo¬ 
tion of a particle is resonant. The treatment of motion is 
2D, and the theory shows that the resonant phenomena 
can be described by a “resonant driving term”, which 
is a quantity that takes into account the distribution of 
the multipolar errors along the ring and characterizes the 
resonance stop-band. 

The interplay of the stability domain and resonances is 
a very difficult subject and there seems to be no general 
theory that explains it completely. However, simulations 
show that in the proximity to a resonance the regime with 
stable motion typically shrinks, and the phase space frag¬ 
ments into regions with very different dynamical proper¬ 
ties: from stable motion, to unstable, and chaotic mo¬ 
tion [T]. 

The difficulty of the dynamics in the vicinity of a res¬ 
onance is remarkable even in regimes in which the mo¬ 
tion do not exhibit chaoticity. Concepts as fix-points 
easily understandable in a 2D phase space assume quite 
challenging aspects when the effect of resonances is con¬ 
sidered in the 4D phase space. For regular motion in 4D 
phase space, the fix-points become exotic objects that we 
call “fix-lines” [8|. The structures created by the 2 de¬ 
gree of freedom resonances have been subject of studies, 
which addressed the problem from the mapping approach 
[9lllQ|. 

Fix-lines become particularly relevant when the non¬ 
linear dynamics is affected by space charge in bunched 
beams. The phenomena of periodic resonance crossing 


induced by synchrotron motion and space charge, create 
new and more complex dynamics which is very signihcant 
for present projects like the SISlOO of the FAIR project 
where bunches are stored for a long term in presence 
of significant space charge m, or for the LIU project 
at GERN, where PSB, PS and SPS are upgraded in in¬ 
tensity and the same topology of high intensity problem 
is encountered [T2|. For ID resonances the interplay of 
space charge with hx-points has been extensively studied, 
numerically and experimentally [Hill]. On the other 
hand, studies on the interplay of space charge and 2D 
resonances has never be attempted because of the re¬ 
markable difficulty in characterizing stable resonances in 
4D phase space (the fix-lines). 

In this paper we investigate the fix-lines in the prox¬ 
imity of a coupled third order resonance and find the 
relation of these objects with the stability properties of 
motion. Hence we extend the theory as presented by 
Guignard [7|- 

The plan of the paper is the following: In Sect. H we 
give an overview of the 2D resonances: this phenomeno¬ 
logical analysis categorizes the resonances and discusses 
the main features of ID and 2D resonances. The concept 
of fix-line is introduced and examples are shown. In Sect. 
HI the phenomenology of the third order coupled reso¬ 
nances is being outlined. Fix-lines and tori cut are dis¬ 
cussed using proper cuts in two dimensional planes. Sec¬ 
tion IV describes the theory of particle motion in a con¬ 
stant focusing channel equipped with one sextupole, the 
theory is developed following a perturbative approach. 
Section V deals with the problem of removing the time 
dependence from the Hamiltonian. We find that there 
are infinite canonical transformations able to remove the 
time dependence. This finding extends the theory of 
Schoch, and in fact it directly leads to a theory of the 
fix-lines (Sect. VI). Section VH treats the issue of the 
stability of the motion around a fix-line, and presents ex¬ 
act formulas for the case of a single resonant term. The 
analysis of stability allows to derive the secondary tunes 
of the motion around the fix-lines. In Sect. VHI we dis¬ 
cuss the stability of motion in the vicinity of the third 
order coupled resonance, and find with large generality 
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a complete characterization of the stability of motion for 
the case of the dynamics dominated by one single reso¬ 
nant term. Tests and comparison with tracking simula¬ 
tions are shown as well. 

In Sect. IX we extend the analysis to an AG structure 
equipped with many thin sextupoles. Following the cor¬ 
responding arguments used for the constant focusing lat¬ 
tice, we directly retrieve a description of the dynamics in 
the proximity of a third order coupled resonance in terms 
of the resonance driving term as obtained from the theory 
of Schoch (for an arbitrary distribution of perturbative 
sextupolar errors). The comparison with multi-particle 
simulations is addressed as well. Lastly, in Section X 
we show that our theory allows to explain the results of 
the tori cut presented in Sect. II. Section XI is devoted 
to the conclusions and Sect. XII holds the acknowledg¬ 
ments. In the appendix A (Sect. XIII) we discuss the 
perturbed solutions around the fix-lines, and in the ap¬ 
pendix B (Sect. XIV) we elaborate in more details on 
mathematical aspects of the discussion in Sect. VIII. 


II. PHENOMENOLOGICAL ANALYSIS OF 
RESONANCE STRUCTURES IN 2 DEGREES OF 
FREEDOM 

A. Categories of resonances in the 4D phase space 

Resonance structures in ID of freedom have been stud¬ 
ied in depth and will not be mentioned except for a 
schematic comparison with the 2D case. 

In 2 degrees of freedom we have basically 3 different 
types of resonances. In Fig. we find an example of a 
vertical ID resonance addressed with “1”. In the vertical 
plane you do expect to see a resonance structure with 
distinct resonance islands, whilst in the horizontal plane 
the motion is non resonant. Therefore we will expect to 
find a number of resonances vertical islands that may be 
widened due to deformations in the horizontal plane. We 
will not follow this type of resonance any further since it 
is basically very similar to the classical and well-known 
ID resonances. 

At the heart of this paper are the coupled resonance of 
type “2” in Fig.j^ The theoretical understanding of such 
coupled resonances including detuning with amplitude 
will be addressed in the following chapters. Here we want 
to show how they affect the motion in phase space. 

A peculiar type of resonance is depicted with “3” in 
Fig. [3 This is the 2D analogon of ID fix-points, i.e. 
these fix-points are equivalent to the central fix-point at 
the origin of the 2D closed orbit if it is different from zero 
in any of the 4 coordinates. To some extent the motion 
in the vicinity of fix-points is better defined then close 
to a fix-line structure. The authors envisage a separate 
report just dedicated to those 2D fix-points. 



FIG. 1: Types of resonances in 4D phase space. Resonance 
lines are typically presented in the Qx-Qy diagram. 


B. Schematic description of 1 and 2 degree of 
freedom resonances 


Figure 13 shows schematically how the 2D phase space 
looks for a stabilized resonance, i.e. the detuning with 
amplitude creates island instead of unbounded motion 
along hyperbolic lines to infinity. It is important to men¬ 
tion that the phase space in the plane of the resonance for 
the type “1” resonance in Fig. [l] would look quite similar 
except that one will experience a widening of the projec¬ 
tion of island motion that in the ID schematic depiction 
is a measureless thin line around the fix-point. 

This classical island structure in ID has mostly regu¬ 
lar motion in the inner part at small amplitude to the 
central closed-orbit, which is the green area in the Fig. [3 
For completeness one has to mention that careful inspec¬ 
tion of the motion in that “green” area will reveal other 
fine island structures. Besides the regular motion around 
stable ID fix-points, depicted as the red area in Fig. [3 
one finds deterministic chaotic motion in the vicinity of 
the unstable fix-points and the separatrix which sepa¬ 
rates the island motion from the “green” area at lower 
amplitude. At larger amplitude outside the island struc¬ 
ture (for simplicity not shown in Fig. [3 one would expect 
regular motion up to a limit where motion becomes un¬ 
bounded. This limit is called the dynamic aperture. 

Notice that there is a dashed circular line that con¬ 
nects the fix-points which is denoted by “undisturbed” 
KAM torus. By that we mean the undisturbed motion 
in absence of the non-linearities of the system. 
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FIG. 2: Features of a ID resonance stabilized through detun¬ 
ing with amplitude. 


The last argument is a good link to resonances in 4D 
phase space, the real emphasis of this work. In 2 degrees 
of freedom the phase space is four dimensional. However, 
the motion of each individual particle is restricted to a 
two-dimensional manifold. In passing, we like to stress 
that the restriction to these two-dimensional manifolds is 
the very reason why there is a fundamental difference be¬ 
tween motion in 2D and 4D phase space respectively: in 
the case of 1 degree of freedom each regular particle mo¬ 
tions is on a curved line that separates phase space into 
an inside and outside part so that no particle can cross¬ 
over from one part to the other. However, in 2 degrees 
of freedom particles move on two-dimensional manifolds 
that cannot separate the 4D phase space into a inside and 
outside part since that would require a three-dimensional 
manifold. In consequence this means that in 4D phase 
space all chaotic regions are linked through-out phase 
space such that there is no longer a strict concept of the 
dynamic aperture in systems larger than 1 degree of free¬ 
dom. 

While particle motion in 2D phase space is pretty much 
self-evident one needs some preparation to finds one’s 
way in 4D phase space. To this end it is helpful to look 
at a 3D projection of linear and uncoupled motion in the 
4D phase space. In Fig. the motion in phase space of 
x,x',y is shown while y' is still hidden. Part a) shows 
x,x',y at some angle around the x' axis: it becomes evi¬ 
dent that the motion is on a torus. When looking straight 
from the horizontal viewpoint (x, x' in part b) one finds 
a perfect ellipse without any thickness and for mixed co¬ 
ordinates (y, x' in part c) one finds a rectangular shape. 
It goes without saying that the ( y, y') phase space is also 
an ellipse and that for all mixed coordinates one finds a 
rectangular shape. 

It is important to realize that we are still lacking the 
4*^ coordinate y'. A good way to describe this extra di¬ 
mension is to understand that the particle moves on a 
double layered torus, i.e. two-dimensional manifolds in 


4D. Introducing non-linearities tends to separate the lay¬ 
ers of the torus so that very complex structures will be 
found in phase space. This becomes evident in Fig. 
which shows a case with strong non-linearities switched 
on: large amplitudes have been chosen but the motion 
remains regular and refined to a two-dimensional mani¬ 
fold. The rotation angle of part a) is chosen to make the 
double layer structure strikingly visible. 


Traditionally, there have been attempts to use other 
means E.G. a color code to allow for a full grasp of the 
dynamics in 4D. Set aside spectacular aesthetic pictures 
effectively this did little to substantially enhance our un¬ 
derstanding of those dynamics. What really counts is 
that one can use cuts in phase space that take away 
overlapping pieces of the 2D manifold in the phase space 
projections. 


Schematically one can have a look at how the “undis¬ 
turbed” 2D KAM torus in case of linear motion is broken 
up into a fix-line structures as shown in Fig.[^ One finds 
that of this torus we are left with a set of alternating 
stable and unstable fix-lines: in the vicinity of the sta¬ 
ble fix-lines the particles move on tori around the fix-line 
(again a 2D manifold), while the motion at the unstable 
fix-lines and the limiting “separatrix torus” is of chaotic 
nature. 


The comparison of Fig. and Fig. shows how island 
structures in 2D are a generalization of those in ID, albeit 
in 4D we have additional complexity. A more complete 
analogon to the ID case is to be expected around 2D 
fix-points as shown below. 


In case of a fix-line there is one condition fixed be¬ 
tween the tunes therefore the stable 2D manifolds break 
up into ID objects. In case there is a second condition 
fulfilled between the tunes these fix-lines are further bro¬ 
ken up so that only fix-points, stable or unstable, remain. 
In Fig. [^it is schematically shown how the particle mo¬ 
tion evolves around stable fix-points. In fact, one expects 
structures that resemble the motion around the central 
fix-point around the closed-orbit, i.e. regular motion in 
small distance to the fix-point and the set-in of chaos at 
larger amplitudes. In the figure the lower left system the 
evolution of the chaotic particles is colored in green for 
better visibility. 
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b.) 



c.) 



FIG. 3: Linear motion in 3D projections of the 4D phase FIG. 4: Non-Linear and regular motion in 3D projections of 
space. The coordinate axis are denoted by 1, 2, 3 which stands the 4D phase space. The coordinate axis are denoted by 1, 2, 
for x,x',y respectively. 3 which stands for x,x',y respectively. 
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FIG. 5: Features of a 2D fix-line resonance. 
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FIG. 6: Features of a 2D fix-point resonance. 


C. Simulation examples of 2D Ax-line and fix-point 
resonances 

In Fig|^and Fig|^we have already shown simulation 
examples for linear and non-linear motion. As explained 
above we are applying cuts in phase space such that one 
can inspect the two-dimensional manifold unobstructed 
by projections of other pieces of this manifold. Fig. 
part a) depicts a regular and dense coverage of such a 
manifold. 

In part b) the amplitudes have been varied in such a 
way that the tunes are on resonance and the motion is 
restricted to a torus around and with a small distance 
to the stable fix-lines. When the distance to the stable 
fix-lines is increased the motion becomes chaotic and the 
unstable fix-lines become visible. In Fig. part c) both 
the stable and unstable fix-lines are indicated with red 
and blue arrows respectively. The two types of fix-lines 
are interleaved and the motion near the unstable fix-lines 
exhibits the typical large chaotic variations. It should be 
mentioned that the motion of part c) encloses that of 
part b). 


a) 



FIG. 7: 2D fix-line phase space cut-out. 


Lastly, we are also presenting an example of particle 
motion in the vicinity of 2D fix-points. In Fig. part a) 
close to the fix-points the motion type is regular. When 
the distance to the stable 2D fix-points is increased, the 
motion will exhibit chaotic behavior which is the case in 
Fig. part b) where one finds fuzzy images due to the 
chaoticity. In part a) the motion around one of the fixed 
point is circled, while the motion in part b) is around 
the same fix-point. It is interesting to demonstrate that 
chaotic motion of part b) encloses the stable motion of 
part a). To this end these two trajectories are drawn on 
top of each other and turned in the three dimensional 
phase space projections x, x, y so that one can see that 
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both motion types evolve around the same fix-point (see 
Fig. § part c) ): one can nicely find regular motion close 
to the fix-point and fuzzy, i.e. chaotic motion at larger 
amplitude. It has to be mentioned that this holds true in 
all 3D projections, albeit without a visual demonstration 
shown here. 


D. Locking to a 2D fix-line strncture 


In the analysis of motion in the vicinity of resonances 
it is important to consider a balance of resonance driving 
terms and the detuning with amplitudes. Often, too sim¬ 
ple models do not really describe what will be found in 
well designed accelerators. Designers go to great pains to 
minimize driving terms and choose tune working points 
with optimal stability, i.e. large dynamic aperture. 

To this end we are using an example with a sizable 
number of FODO structures each with sextupoles with 
an average non-zero value and some random component. 
This set-up can be seen as a reasonable good approxima¬ 
tion for a realistic accelerator. 

To study the effect of the sextupole coupling resonance 
Qx + 2Qy the bare tunes are set close to that resonance. 
We have tested II runs incremented by a fixed 4D coor¬ 
dinate vector. 

In Fig. Impart a) the phase advance per turn is shown 
over several thousand turns. Despite the apparent large 
variations of the phase advance one also finds the tunes 
of each case and a line connecting the tunes. Due to 
the specific detuning with amplitude the resonance line 
Qx+‘2Qy (shown in the figure as well) is being approached 
until at step “8” when the motion is locked to the res¬ 
onance. With larger amplitudes in this 2D case we are 
finding a phenomenon not known in ID: for the next 3 
steps the tunes are continuing to move but remaining 
locked to the resonance. It is important to understand 
that at all amplitudes in between case “8” and “11” the 
motions remains locked to this resonance. Each point 
though on the resonance line is a fix-line system in it¬ 
self, as described in the previous sections: a fix-line at 
the center and motion around it while keeping identical 
horizontal and vertical tunes. In this study we have en¬ 
countered rather accidentally particle “10” which fulfills 
almost perfectly the condition for a fix-line. 

Figurej^part b) shows a blow-up of this fix-line motion 
in the plane of the phases. It might be surprising at first 
sight that the particle is never actually on the resonance. 
The averaged phase advance, i.e. the tunes are indicated 
by the blue star in the figure. The tunes clearly never 
coincides with the individual phase advance values per 
turn (shown as red dots). 


a) 



b) 

2.0 


I .□ 


0.0 


-I .□ 



c) 



FIG. 8: 2D motion in the vicinity of a stable 2D fix-point. 
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b) 



FIG. 9: Phase-advance per turn; Tnnes (= averaged phase 
advance); Orbiting aronnd the 2D hx-line from case 8 through 
case 11. Part b) is a blow-up of case 10 that is an example of 
a fix-line in the plane of the phases. Interesting enough the 
phases never really respect the resonance condition instead 
the tunes do. 


III. LINKING THE PHENOMENOLOGICAL 
PHASE SPACE DEPICTION WITH THE 
THEORETICAL DERIVATION 


The theory of fix-lines resonance structures will be de¬ 
rived in the following chapters. We would like to use 
a particular theoretical prediction to demonstrate how 
the particle motion of a large scale + 2 x fix-line 
structure due to a single sextupole in a constant focusing 
lattice can be viewed in the 4D phase space. 

It will be shown that in the mixed phase space plane 
y' versus x' we expect a figure 8 like structure (see 
Fig. 10) for the particle motion on the fix-line. 

The particle motion on the stable fix-line and on the 
tori around it creates a fix-line structure that is distinct 
from regular motion around the closed orbit. To visualize 
this fix-line structure we show the stable fix-line and 3 
tori around it in a restricted phase space as described 
in the earlier chapters. We have to do so to view only a 


single passage of the fix-line structure in the y' versus x' 
phase space projection. To this end we restrict x and y 
to positive values. Figure 0 gives a first impression of 
this fix-line structure. However, it remains difficult to 
figure out graphically how the tori are stacked around 
the fix-line. Therefore, we apply a more vigorous cut in 
the phase space by restricting the particle motion to a 
very small angle in either the horizontal or vertical plane 
and depicting the resulting phase space projection in the 
other plane. 



FIG. 10: The fix-line is seen in the y' versus x' plane as a 
figure 8 like structure. 



x’ X 10 


FIG. 11: Single passage of the fix-line structure in the y' 
versus x' phase space projection. 


Figure shows that in the horizontal plane the fix¬ 
line (red dots) lies indeed in the center with the particle 











motion located on the three tori with increasing distance 
around the fix-line. The same is true in the vertical plane 
as seen in Fig. |13[ except now we find two such struc¬ 
tures which reflects the fact that we are investigating the 
Qx + 2 X Qy resonance. 



X y. 


FIG. 12: Horizontal phase space projection after restricting 
the vertical motion to a small angle. 



y xlO 


FIG. 13: Vertical phase space projection after restricting the 
horizontal motion to a small angle. 



FIG. 14: Right part (see Fig. 131 of vertical phase space pro¬ 
jection after restricting the horizontal motion to a small angle. 


We have added a close-up (Fig. 14) for positive y values 
so that details with the tori and the fix-line at the center 
become more apparent. 

In the following chapters we derive a theory of the par¬ 
ticle motion in the vicinity of the third order coupling res¬ 
onance. We start from a constant focusing model, which 
allows a simpler description of the dynamics. Starting 
from the description of the linear motion, we derive the 
motion of a particle by varying the constants (emittance 
ttx^ay and phases ipx,^Py). We obtain a system of dif¬ 
ferential equations in ax-,ay,ipx-,fy which is “time de¬ 
pendent”. Then we demonstrate that there are infinite 
canonical transformations capable of removing the time 
dependence in the Hamiltonian of the system. Once the 
time dependence is removed, the simplest solution is ob¬ 
tained when the new variables are constant. This par¬ 
ticular solution is proven to be the fix-line, and we will 
obtain analytically the result of Fig. 10 It turns out 
that the parametrization of the canonical transformation 
is the parametrization of the fix-lines, hence we find that 
there are infinite fix-lines contrary to the ID case with 
one set of only 3 (unstable) fix-points. Perturbations 
from the solutions allow to discuss the stability of the 
fix-line and to derive secondary tunes, i.e. counting the 
oscillations around the fix-lines. Moreover, we find that 
sextupoles by themselves are capable to create stable fix¬ 
lines, a feature not know in the ID case. 

The particle dynamics become better described in 
a special combination of the “varying constants” 
axjCLyjipxT^Py- we find that the motion in this special 
combination has another invariant that allows to discuss 
the stability. Consequently, we are able to describe the 
stability domain in the vicinity of the third order cou¬ 
pled resonance. As may have been expected we find that 
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the unstable fix-lines define the border of the stability 
regime. We derive the properties of the stability domain 
in terms of scaled quantities. The scaling factor depends 
on the distance from the resonance and on the driving 
term. 

Lastly, we apply our approach to the AG machines. 
The conclusion remains the same, and the scaling pa¬ 
rameters is incorporating the complexity of the AG lat¬ 
tice along with the distribution of sextupolar errors. 


IV. EQUATION OF MOTION FOR LATTICE 
WITH CONSTANT FOCUSING 


Let’s start with the equations of motion describing the 
dynamics of a particle in a circular accelerator of radius 

R. 


^ +k^x = foo{s,x,y), 
^ + fcyj/ = fy{s,x,y), 


( 1 ) 


and 


H. = J - ^) ■ (6) 

Therefore for a localized sextupole of integrated strength 
K 2 located at d = 0 in a circular machine, using Eq. 
we write 

k2{0)= -^d(6»-b27rm). (7) 

m—— C50 

where the sum over infinite delta functions (5() is used for 
modeling the periodicity of one sextupolar error in the 

machine. The Hamiltonian Hi becomes then 

Hi=RK2 Y + 

m— — oo ^ ' 

The solution of the unperturbed system with Hamilto¬ 
nian Hq is the solution of the canonical equations 


here kx,ky are constants. The tunes of the machine 
QxiQy are defined as \/k^R = Qx, \/T^R = Qy In 
our discussion fx,fy are created by a magnetic multi¬ 
pole, then 

- fx{s) + ify{s)=^^^^{x + iyY. ( 2 ) 

n! 


dx 

dHo 

dx' 

dHo 

d9 

dx' ’ 

d9 

dx 

dy 

dHo 

dy' 

dHo 

d9 

dy' ’ 

d9 

dy 


which yields the equations of motion, that are 


With this notation for n = 1 we get —/x(s) = 
ki{s)x, fy{s) = ki{s)y, hence, a positive fci yields a 
focusing force in the x plane. 

Let’s make the change of the variable s to the new 
variable s —?► 0 so that ^ We use the notation 

x' = ^. With this change of variable the equations of 
motion become 

x"+ Qlx = R‘^fxis{9),x,y), 

y” + Qlv = R^fy{s{d),x,y), 

and the length of the one turn is A0 = 27r. The inte¬ 
grated strength of ^ 2 ( 0 ) in the 9 coordinate is linked to 
the integrated strength in the s coordinate via the rela¬ 
tion 

/ k2{9)d9= k2{s)^ds=^ (4) 

dSi Js^ as R 

where K 2 is the integrated strength in the s coordinate. 


A. Perturbative approach: the dynamics with one 
sextupole 


We will consider next the accelerator with a single sex¬ 
tupole for exciting sextupolar resonances. The Hamilto¬ 
nian of this system is H = Hq + Hi with 


rr I '2 , Qx 2,1 '2 

ffo = 


Ql 2 

— V 
2 ^ 


x" + QxX = 0 , y" + Qly = 0 . ( 10 ) 

The solutions are 

X — \/jdxax cosi^Qx^ T y^x')i (11) 

y — ^l^yay cos(^Qy9 ipy'j^ 

with /3x = 1/Qx,l3y = ^IQy The parameters ax,ay are 
the single particle invariant, i.e. the usual emittances ex¬ 
pressed in the reference frame with 9 coordinate. Hence 
they are not the emittances in the laboratory frame. 

Next we use the canonical theory of perturbation, in 
which we add Hi and compute the evolution in the sys¬ 
tem Hq -|- Hi of the otherwise constants ax,^Px,ay,ipy. 
The equations of motion of ax,<Px,ay, (py are 


a 


I 

X 


a 


/ 

y 


dHi , JHi 

dpx ’ dax ’ 

dHi , dHi 

dipy ’ day 


( 12 ) 


Note that this is a canonical system with Hamiltonian 
2Hi. Next we substitute in Hi the Eqs. E] re-written in 
the form 


— \/%x 


y = \/% 






( 13 ) 


(5) 


2 
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with Ax = QxO + '^x, and Ay = QyO + ipy defined just 
for convenience. We also expand the periodic S as 

OO -| oo 

S(0 + 27rm) = — (14) 

m— — oo 1— — 00 

and substitute it in Hi. After some algebraic calculation, 
Hi becomes 


V. REMOVING THE TIME DEPENDENCY 


We redefine the variables ipx , ifiy so that they incorpo¬ 
rate the term Ar9. Let’s define 


— ^x ix9: 

‘fy ~ ‘Py 


( 22 ) 


Hi =RK2^{l3xaxf^'^ X 

3 oo .V 

H ( .) cos[((2j-3)g; 

j=0l=-ao 


1)6 + (2j - i)^x\- 


l I _ 00^1 

~ IOtt 2 ^I X! (./) 

[/=-oo j=0 j'=0 ' 

cos[[(2j^ — l)Qx + (2j — 2 )Qy — l)] 6 + 

+ {2j — l)'~Px + {2j — 2)ipyW . 



(15) 


At this point we keep the slow varying terms only, that is 
the terms close to the resonance Qx -I- 2gj, = N. Here N 
is an integer used throughout this article with the same 
meaning: it defines the location of the resonance. Let’s 
define 


Ax — —N -\- Qx -l- 2Qy, (16) 


with tx,ty some constants to be defined later. Then the 
argument of the cosine in Eq. becomes 

Ax6 + (fix ~\~ ‘2-^Py = (px ‘2-<Py + [Aj. — [tx + 2ty)\9j (23) 

therefore if A,, — {tx + 2ty) = 0 the time dependence dis¬ 
appears. We then consider the coordinate transformation 

5 

Px = Px + tx9 ^ (24) 

Py ~ Py ^y^- 

If we define 

Rsi{Axi Pxi ^yi Py^ 9) — Hgi (Ua,, Px i'x9^ 0,y^ Py ty9j 9) , 

(25) 

then 


then the slower varying terms are obtained in the second 
term for I = N,j' = l,j = 2, and I = —N,j' = 0,j = 
0. Hence the slower harmonics Hamiltonian Lfsi, which 
mainly contributes to the dynamics is 

Hgi — yy 'x/ Px^xi^y^y cos )Ax9 -\- x 2p>y). (17) 

iOTT 

We define for convenience 

^ = -RK2^^xPy. (18) 

If instead we consider a continuous sextupolar component 
with strength k2(9) = JC2COs{N9), the expression for A 
is 

A = (19) 


2?A1 = 2?^ = -< = -i’,. 


dpx 


dpx 
.dHsi 


dHgi laaigi , 

^ dax ~ dax 


,dHsi 




2—= 0,y — Qiy, 


dpy dpy 

,dHsi x,dHsi 


da,. 


= 2 - 


da,, 


= p'v = p'y - h- 


(26) 


Therefore we have constructed a function (Eq. 25) that 
is time independent, but if we try to derive the equations 
of motion from it, we find Eqs. |26| These equations are 
clearly not the canonical equations, that means that the 
function Eq. is not an Hamiltonian. This happens 
because in 


Note that JC 2 in this experssion is the strength of a dis¬ 
tributed sxtupole not to be confused with the integrated 
strength K 2 in Eq. We therefore find that the general 
form of the slowly varying Hamiltonian is 

Hsi = A^/aQuy cos[Ard Px + “^Py] (20) 

and that it is time dependent as 
OH 

—^ =-A^/axaySin[Ar9 + Px+ 2ipy]Ax ^ 0 . (21) 


= (27) 

there is the constant term —tx that should not be there. 
We then define the new function 

R-Slipxj Px, dy, Py^ 9) =Hsl {jXx, Px ix9, (ly, Py ty9,9) -)- 

_ tx _ tx 

+ ax— + ay — 


(28) 
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which is also time independent. Next we check if one can 
get the canonical equations 

5^. dip. 

dH,i dHsi , 

- \- tx — Px + tx — Px, 


da. 


da. 


dHsi dHsi , 

2 _= 2 _= —h 

dpy dpy 


2^ = 2 

da„ 


y 

dHsi 

day 


+ iy — P'v — ‘^'v 


(29) 


Therefore we obtain 


= ‘fx 


(30) 


dpx 
^ dHs i 
ddx 

2^ - -a' 
dpy - 

^ ddy 

which are the canonical equations with the Hamiltonian 
2Hsi. By construction the Hamiltonian Hgi is time in¬ 
dependent. In fact 


Hslidx,Px,dy,Py,0) = 


dxdy COS[<Pa; -I- 2py\P 

tx _ (31) 


a. 2 


a 

ay 2 


does not depend on 9. Note that the procedure de¬ 
scribed here is equ ival ent to transform the Hamilto¬ 
nian 2Hsi (from Eq. 20) to the Hamiltonian 2Hsi (from 
Eq. 31) through the canonical transformation Eqs. 
constructed by using the generating function 

^2^Px: Pyi ^x: ^y) — Px^x “t“ Py^y T tx^x9 -)- tyQyO (32) 


with tx -l- 2ty = Ar- 


A. Summary 


We can summarize the situation as follows. There are 
infinite canonical transformations of the form 


Px - Px ^x9 ^ 

~ ‘fy ~ ^y^j 


(33) 


xy, 


with tx + 2ty = Ar that allow to describe the time evo¬ 
lution as solution of canonical equations of the time in¬ 
dependent Hamiltonian 


Hsi (a.: Pxi ; V^y) — ^x^y 0Os[(/^3^ -)- 2 (pyj-)- 

~ tx , ty _ (34) 

CL. — + fly TT + -^(fla:, fly)- 


The function J^{dx,dy) collects all the nonlinear terms 
arising from the theory when other multipolar compo¬ 
nents (non resonant) are included in the Hamiltonian Hi. 
Hollowing the procedure previously outlined, the slowly 
varying terms will be found in the zero-th order harmon¬ 
ics, that will create an additional term in the slowly vary¬ 
ing Hamilton Hgi- Space charge can be treated in the 
same way. The explicit form of the canonical equations 
in the coordinates dx,Px,dy, py is 

-d'x = -2Ai/^dy sin[(^a; -|- 2py] 

2v«x (35) 

-d'y = -AAi/^dy sin[(^a; -|- 2py] 

Py — 2A\Pax c.os\^x T ‘^Py\ ty ~\~ 2jFy 

where we used the short notation J^x = dTjddx-, ^^nd 
Ty = dTjddy. Erom these equations we find immedi¬ 
ately with 2d!x — d'y = Q : 

2dx — dy = C (36) 

with C a constant. As the constants tx,ty must satisfy 
the condition tx + 2ty = A., the solution dx, p.^dy, Py 
depends on two free parameters, a choice could be t., C. 


VI. SELECTING THE CANONICAL 
TRANSFORMATION 

In this section we search in the space (d., Px,dy, Py) 
the solution of = 0,Px = 0, d'y = t),p'y = 0. The 
canonical equations then read 


0 — Q ^ — A-^ a.ay sin(c^a, -)- 2py^^ 

^ dHsi . 1 _ r ~ - 1 tx .r- 

0 = = A—^ay cos[(/ 73, -h 2py\ + — +Fx, 

OQjx ^ V ^x ^ 

0 = „ f = -2Ai/^dy sin[(pa; -|- 2py], 


0 = 


dpy 

dHsi 


= Aa/^ COS[(^;j -I- 2py\ + ^ +J^y. 


ddy "V > -^yj W 2 W.y. 


(37) 


The equations [37|l|37|-3 are satisfied for 

sin[(^2, -I- 2py\ = 0, (38) 


that is for px + 2py = ttM, with M = 0,1. The other 2 
equations become 


1 


0 = A^^d^(-l)'^ + ^+.F„ 

/ax 

0 = A-\/^(—1)^^ + "^ + -^y 


(39) 
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From the condition tx + 2ty = we can eliminate tx 
and we get the system 

f 0 = A ay[-l)’^ + - ty + 

J 2y'ax 2 

0 = Aa/^(-1)^ + Y + 

This system can be solved in ax,ay, from which we find 
also C through Eq. We have therefore established the 
correspondence 

ty y (da; , Uy , fa; , d/) , 

which means that all the solutions of Eqs. are on 
a 1-dimensional curve in the plane (da;,dy). This 1- 
dimensional curve can be expressed as a function of 
dx,ay, without any reference to tx,ty,C, and together 
with the phase relation one gets 


\ 0 = A(-l)“ 

dy 

_L 9, 


\ ^x 

( (^3, 4- 2(fy = ttM. 



Ar 

+ ^+^x 




(42) 

The significance of this equation is the following: for 
any pair (dx,ay) that satisfies Eq. 42 (top) there exists 
an associated canonical transformation characterized by 
the pair {tx,ty) obtained from Eqs. |39[ In the new coor¬ 
dinates identified by {tx,ty) the pair^a;,dy) joined with 
(fx + 2(py = ttM solves Eq. i.e. in the system of 
coordinates identified by (tx,t^, {dx,dy) are constant. 
Therefore all the {dx,dy) satisfying Eq. H^are all possi¬ 
ble stationary solutions of Eq. 37 (for the correspondent 
canonical transformations). The set of these solutions 
(dx, dy) is a 1 dimensional curve (or a collection of 1 di¬ 
mensional curves for complicated Tx,J^y dependencies). 


A. The fix-line 


Next we discuss the meaning of the solution d^ = 
0, ‘fx ^0,d'y = 0, f'y = 0. As f'x; = !f'y^0we get 


— 0y,Oj 


(43) 


and returning back to the coordinates of the Hamiltonian 
Hsi (Eq. 20 there written without T) we find 


‘fx — ‘fx,0 ^x^: 

fy ~ fy,o ~ tyS, 


(44) 


with ‘fxfi + 2:^y_o = ttM, resulting for the coordinates of 
Eqs. [TT1 


^ — \/ Px^x COs[(d^a; ix')^ 4“ ; 

^ — i^Qx tx^ \/Px^x Sin[(d^a; ix^^ 4“ ^a;,o]; 

y \/Py0>y C0s[((5y 4- 9^y,o]; 

y (.Qy ^y) \/^y^y ^y)^ 4- ^y,o]' 


(45) 


We will next prove that the points of this trajectory at a 
given longitudinal position (for example at 0 = 0) lie on 
a one dimensional closed curve that we call fix-line. The 
proof of this follows: 

By using the conditions ifx+2ipy = ttM, tx+2ty = Ar, 
and Eq. [T^ it is straightforward to find that 

{Qx—tx)0 + ‘fx,O + 2[{Qy—ty)0 + fyfi] = N O + TT M. (46) 


The Poincare surface of section is identified by the con¬ 
dition 9 = 2TrAf with JV an integer correspondent to the 
A7-th turn. Next we limit the discussion to the a:, y co¬ 
ordinates without losing generality. By using Eq. we 
find that on the Poincare surface of section the coordi¬ 
nates X, y at the turn M become 


Xj\f — \/ i^x^x COs[((5a; tx)27rAf 4- 

yjf C0s[(f^y ty'^2'ITj\f 


and using Eq. 46 we find 

xj\f — \/Px^x cos[ 2((^y ty)27rA7 
— 2(^y Q Tribf -f A27rA7], 

yjf — 's/i^y^y COs[(f^y ty^2lTj\f V^y,o]; 


(47) 


(48) 


with <fy^o Ets an initial phase that sets the starting point. 
Therefore we reach the result that the coordinates x, y 
of Eq 1^ on a surface of the Poincare section, can be 
parameterized (ignoring multiples of 27r) as 

xt = yjPxttx cos(-2t -f ttM), 

- (49) 

Vt = V/Syfly cos(t), 


a) b) 



. -75 .. 

-4 -2 0 2 4 , -4 -2 0 2 4, 

X X 10 x' X 10 


FIG. 15: The fix-line as predicted by the theory (red dots) 
and by simulations (black dots). Pictures a), b) circular shape 
shows that dx,dy are constant. 
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where Qx, ay must satisfy Eq. 42 (top), and 0 < t < 27r 


is here a variable that parameterizes the curve. Therefore 
the coordinates x,y of the solution in Eq. |^are always 
found on the curve expressed by Eq. 1^ The same argu¬ 


ment applies to any other pairs of coordinates of Eq 45 


We therefore find that at a given longitudinal position 
the trajectory Eq lies on a closed curve. This curve 
is invariant in its shape, and each point in this curve 
is mapped after one turn again on the same curve, for 
this reason this curve is called fix-line, which is a gen¬ 
eralization of the fix-point for ID stable resonances. In 
Fig. [^we compare the fix-line as predicted by the the¬ 
ory (red dots) with simulation results (black dots). The 
figure shows 4 pictures for the planes, x — x',y— y',x— y, 
and x' — y'. It is clearly visible that for the variables 
x, y, and x' — y' the shape is that one predicted from 
the theory. However, the planes x — x'^ V ~ v' exhibits 


a circular shape showing that dx,ay are constant! The 
simulations are obtained by exciting a single harmon¬ 
ics, and Qx = 1.324, Qy = 1.84, = 4.64 x 10“^, and 

A = -7.444 X 10 


-3 


'y 

,-i 


B. Stability 

We consider here a particle slightly displaced from the 
stationary solution of Eqs. 37 (this solution is the fix¬ 
line) of the amount Sdx, Sdy, S(px, 5(py. In order to obtain 
the equation of motion of the perturbation we expand 
equations 35 in Sdx,Sdy,6(pxj^‘Py, and keep the first or¬ 
der only. The canonical equations of the perturbation 
take the form (properly re-arranged) 


Sd' = -2SA^5dx - 2-^^^"^ 


difxddx 


difxddi 


Sdy - 2^^-S^S<fx - 6<py, 


d(pl 


d(pxd(py 


^ a ~ ^ a ~ O'- ^ o - o - 2 

^^y^^X UiPyUQjy U<^yU(Px 


5(p'^ = 2^-£^6dx + 2 Sd,. 


<5^; = 2 


ddl 

d^Hsi 

ddyddx 


ddxdifx 


Sdx 


ddxddy 

2^^£^5dy + 2.^r^^5fpx 


2 ^^ .... diPx + 2 ^. diPy, 


ddl 


ddydipx 


ddxdfpy 

ddydify 


5^y. 


(50) 


r 


where the slowly varying Hamiltonian is given by 
Eq. |34| Any term that has only one derivative with 
respect to (px , fpy yields zero because that term is pro¬ 
portional to sin[(^a; -|- 2(py\, which is zero on the fix-line. 
Therefore we get 

5d'^ = -2ZS^Sp>x - 2 ^!^"! Sify, 




d(pxd(py 


oa„ = -2 „, difx - 2 dify, 


dfpydfpx 


Sdx 


ddl 

d^Hsi 

ddyddx 


We define for convenience 


Sfp'x = 2 


S^'=2 


Sax 


difi 

9 A' 

ddxddy^ 


A = 2Si=2l-"^'^- 




2 d(pxd(p 


4 dp=l ■ 


Axx — 2 


d^Hsi 
ddl 

A -2AAl 

■^yx — O'- o~ 

oa.iOar, 


A -2^ 

■^xy — ^ o" 

oa^oa. 


•y 


A — 9 

J-Lyy — ^ 


d^Hs, 

ddl 


(51) 


(52) 


(53) 


Therefore we always find that the motion of the pertur¬ 
bation is given by the following system of equations 

Sax — S^Sp?x 2\Sp^yj 
Sd'y = —2\S(px — ‘i\S(py, 

S(Px — A^XxSdx -^xySdy ^ 

SiPy = A^yxSdx ^yySay 

which after some algebra becomes 

Sdx ~ ^ Sdx 4“ A(Aa,y -t- 2Ayy^Cp^ 

Sdy = — up' Sdy — X{AxX + ‘2:Ayx)Cp, 

SSp"x = ~ X{Axx + 2,Axy){S(px + 2Sipy), 

Sfp'l = — X(Axy 2Ayy)(Sipx 2,S(py)^ 
where we defined 

UP = Xi^AxX + 2,Ayx 2Axy -|- 4Ayy). 


(54) 


(55) 


The coefficient Cp is a constant obtained from integrating 
Eqs. 53 we call it Cp with index p as we refer to the 


perturbation around the hx-line. 

The value of Cp depends on the initial condition of the 
perturbation 5dx,Sdy,S(px,S(py. From Eq. 53 we find 

2i5aa; = Sdy Cp. 


(56) 
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We also find 

+ ‘25(py)" =-U!^{6ipx+ 25(py). (57) 

We conclude that for 

> 0 (58) 

the motion of Sax,Sdy, and + 2d(^j, is stable. The 
evolution of 5^px,5(py instead is not bounded if Cp 7 ^ 0 
as a term linear in 9 appears in the solution of the 
equation of motion. The derivation of the solution is 
reported in Appendix A. An interesting consequence 
in the case of Cp 7 ^ 0 is that there exists a neigh¬ 
bouring fix-line identified by d%,dy, and a point on it 
d%,dy,ip%,ify, so that the initial perturbation assumes 
values ((5da;)i, (ddy)i, (^<Px)i, iS(fy)i, and = 2{Sdx)i - 
{6dy)i = 0. That means that the initial perturbation os¬ 
cillates now around the new fix-line point. The proof of 
this result is shown in the Appendix A. 


with 0 < r < 1 and the solutions for the fix-lines take 
the form 



A2 

r 

16A2 

A? 


(1-r)^ 


4A2 


r(l-T) 


(65) 


At this point, as a last step, it should be checked if the so¬ 
lution really exists. Substituting the previous equations 
into Eq. we find that the solution exists as long as 


A Ar 


(-l)^-f-P 1 = 0, 


( 66 ) 


which can be satished only if AAr(—1)^ < 0. Our anal¬ 
ysis shows that in proximity of the coupled 3rd order 
resonance there exist an infinite number of hx-lines. 


VII. EQUATION OF MOTION WITH THE 
RESONANT TERM ONLY 

We discuss the fix-lines for the simplified case of = 
0 , i.e. in absence of detuning induced by other non- 
linearities like space charge, or strong nonlinear compo¬ 
nents due to magnets. The equation of stationary solu¬ 
tions reads (see Eq. 40 with = J^y = 0) 


0 = Av^(-l)^ + |. 


(59) 


We remind here for convenience that A^ = Qx + ‘2Qy. 
Adding the second equation multiplied by 2 to the first 


one ty can be eliminated and we obtain 


0 = A^A(-l) 


M 


A2 

+ ^- (60) 


L2v^ " 

Therefore if the solutions dx, dy exists, it also holds 

ArA(-l)^ < 0 , (61) 

which puts a constrains on tx , ty as follows 


0 = A,A;^dp(-l)“ + A,| 


0 = A, 


d,(-l)^ + A.| 


Afta; > 0 
Aftp > 0 


(62) 

As a consequence if the solution exists, then it must be 
A^ > l^Ax > 0. (63) 

In that case we can parameterize {tx,ty) as follows 

tx=TAr, ty = ^{l-T), (64) 


A. Analysis of the stability around the fix-line 

By using Eqs. the stability coefficients of Eqs. [52| 
are 


- 2A , CL 


l--3/2~ / tx 


,(-!)“ = 


‘2 q,x 


— 1/2/ i'\-^ _ 


Axy = Ayx = 2A-dx ' (- 1 ) 
Ayy = 0, 

A — 2A\J (XxCLy ( 1) — ^yty 


(67) 


where we have used the relation Eqs. |37| on the fix-line 
for = 0. After some algebra we find uj^ = — 'itxty, 

which, in terms of the parametrization of the fix-line, 
becomes 


U) 


= AIt{3t-2). 


( 68 ) 


This equation shows that with r > 2/3, is positive, 


hence Eqs. 54 yields stable motion, whereas for r < 2/3, 


is negative which means that the motion is unstable. 
Therefore this analysis shows that with r > 2/3 a fix¬ 
line is stable. 


The situation is summarized in Fig. 16 top, where 
the collection of all fix-lines is represented. The point 
of merging of the stable with the unstable fix-line has 
= 0 and that corresponds to t = 2/3. 

Figure bottom shows the comparison of the oscil¬ 
lation of an initial condition slightly off a fix-line. It is 
evident that the center of oscillation is not ( 0 , 0 ) because 
2ddx 7 ^ Sdy (hence Cp 7 ^ 0). The red curve is drawn 
from the theoretical model by fitting the amplitude, how¬ 
ever, the wavelength is given by w. The parameters 
of the simulations are the same as for Fig. namely, 
A,. = 4.64 X 10"/ A = -7.444 x lO'^ m"!. The fix-line 
is taken for r = 0.7, which implies tx = 3.248 x 10“^, 


and ty = 6.960 x 10 


-4 
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a) 




FIG. 16: On the top are shown the set of stable and unstable 
fix-lines. The bottom picture shows the oscillations of 
around the fix-line value. In this example Cp 7 ^ 0, hence the 
oscillation of the perturbation Sd^, is not centered around 
zero, i.e. on the unperturbed fix-line. The center of oscillation 
Sd^ is located in this example at —0.2 x 10“®. The black 
curve is obtained from computer simulation of the accelerator 
model, while the red curve is obtained from the theory. The 
wavelength of the red curve is obtained from Eq. | 68 | 


VIII. STABILITY OF MOTION 


We make use of the properties of the equations of mo¬ 
tion Eqs . 35 In particular of the invariant 2dx — dy = C 
(Eq. 36). By using it, we can drop one of the two equa¬ 


tions of Eqs. and summing the other two we get 

— 2i\^ Qjx{2 Qjx G) sin(g?a: -)- 2:^yj, 

(p'^+2ip'y= 2k — ^-={2dx-C)cos[(p^ + 2(py] + 


+ 4AV dx COS[:Pa; -f 2(py\ + Ar -f 2^ j, -)- y. 

(69) 


These equations show that the natural variable for the 
system is 17 = (^a, -I- 2(py, and the equations become 


-d'y. = — 2Ay/d^{2dy: — C) sin[0] 
1 


O' =A 


:(2da, ~ C) -|- 


.yo-x 
A.r 2J-X 


cos[0]- 


(70) 


-My. 


Eor the case in which we take only the resonant term we 
have Tx = 7-y = 0, and the equations take the form 


d!y. = 2KyJ~^(2dx — C) sin[0], 
1 


O' = A 


--{2dx — C) + 4 -\/^ 


cos[0] -b A^. 


(71) 


Now consider the function 

/(dx, O) = 2K\/dx{2dx - C) cos[0] -|- d^Ar, (72) 
and note that 


=A 


1 


:(2dx ~ C) + 


dl 
ddx 
dl 

— = - 2K^/^{2dx - C) sin[0] = -d'y.. 


Therefore we find 


COs[0] -b Ar = O', 


(73) 




That means that /(dx,0) is an invariant of motion and 
the trajectories of particles are set by the value of the 
invariant (that we call Iq). By using /(dx, O) it is possible 
to discuss the stability of the motion of particles. In fact 
given the initial conditions of dx,dy,0 we find C from 
the relation 2dx — dy = C, hence we find the value of the 
invariant /q = /(dx,0). The trajectory of that particle 
is then given by the set T/^ = {(dx,0) : /(dx,0) = /q}. 
We observe that if the trajectory extends up to infinity, 
it means that dx —t 00, hence dj, —>■ cx), while 0 < O < 
27r. In this situation the first term of Eq. [7^ diverges as 
4Adx^^ cos[0], and this divergence is faster than the one 
of the second term, which diverges as dxA^. That means 
that n must change accordingly for keeping unchanged 
the value of the invariant, therefore for dx —t 00 we find 
17 —>■ 7r/2 -b TT-TT, with n an integer. 

We observe that the function / = /(dx,f7) defines a 
function dx = dx(f7), which has the following property: 
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if ri = 27rn (n is an integer) and (ix is not a fix-line, then 
= 0. In fact, as / is an invariant, dI/(Kl = 0, but 


dl dl ddx dl 
dfl ddx dn dfl 


(75) 


that is 


dl _ ddx 


= 0 . 


(76) 


Therefore for = irn it holds sin[n] = 0, hence from 
Eqs. we find d(, = 0. On the other hand, if dx, O is 
not a fix-line then Q' ^ 0 (because (fx = <fy = 0is the 
condition for finding a fix-line). Therefore if O = 27rn, 
and dx is not a fix-line then 


ddx 

dfl 


(77) 


This property tells us that the implicit functions dx = 
dxi^) defined by 7 = I(dx,^) are periodic in O and at 
the periodicity they merge smoothly (that is the first 
order derivative in O is zero at O = 27rn). 

As the fix-lines are described by Eqs. we will con¬ 
sider the following scaling 


- _ A2 ^ 

- _ A2 , 


(78) 


The relation 2da, = dy + C becomes 

^ 16A2 

2Cl£c - Ojy O. 


(79) 


It is convenient to define the parameter 

16A2 C 




/S.I 2 


(80) 


so that we obtain the relation 2dx — dy = 2^. The in¬ 
variant I defined in Eq. rescaled in the “normalized 
coordinates” dx, dy reads 


I — — /r dx ijdx ^) cos n -t- dx , (^^) 


with 


A* = 


A |A,| 
|A| A, • 


(82) 


This relation satisfies the relation /i^ = 1. Note that, 
as 2dx — dy = 2^ it follows that dx must simultaneously 
satisfy the following conditions 


( dx>0 

\ dx ^ C 


( 83 ) 


which express simply the condition da, > 0, dy > 0 intrin¬ 
sic in the definition of dx,dy. In addition we observe that 
all the possible curves I{dx, 71) = const, never intercept. 

The trajectories can be found from 72 as function of dx 
inverting the Eq. obtaining 


cos [71] = fj,- 


1 I dx 

f^x dx - C 


(84) 


where we used the property fjr = 1 (from Eq. 82). To 
study this equation we make use of an auxiliary function 
defined as 


i^(da;) = —AiCOs[72] =-(85) 

ydx C.X s 

This function is defined in the plane d^,, 1^, which is the 
plane about which we will discuss next. Note that ^ is 
here a constant of motion defined by the initial condi¬ 
tions. Hence by varying I we identify a specific level line 
of the invariant. If —1 < < I then ‘3^ = —^cos[72] 

can be solved and we find 72;if^^ <—lorI<^ then 
3^ = —/iCOs[72] cannot be solved so that these values of 
3^ are not acceptable. 


A. Fix-lines and invariants 


Next we prove the following property: consider a level 
line identified by 7, of the form 3^ = 3^{dx,I,0 defined 
by Eq. 85 If this curve is tangent to the line 3^ = I (in 


the plane dx, 3^) then the tangent point {dx, 3/) identifies 
a fix-line. 

Proof: Firstly, we observe that Eq. is simply a 
rewrite of the Hamiltonian, from which we can obtain 
the equations of motion as 


A 


dl 

n' = Ax — . (86) 

dux 


The function 3^ = 3^(dx) is an implicit function of the 
form 7 = I{dx,3^{dx),C} with 7 = const.. Therefore 

d f,. .o... X .X dd dl d3/' ^ 

^d{ax,3^{ax),0 --0. (87) 

If for dx the level line 3^ = 3^{dx) is tangent to a hori¬ 
zontal line, then we have d3^{dx)/ddx = 0. This applied 
in Eq. implies that on the tangent point one finds 


dl 

ddx 


{dx,3^{dx),0 = 0 , 


( 88 ) 


which via the second part of the Eqs. [^yields 72' = 0 (or 
3/'' = 0). Secondly, we use the condition that 3^ = 3^{dx) 
is tangent to = 1; that means that cos[72] = —fi, 
i.e. sin[72] = 0, which yields dl/dVL = 0. Therefore 
via the first part of Eqs. 86 we find d(, = 0. If at any 


time d^ = 71' = 0 then it follows that dx = const., 72 = 
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const., because d,.,!! are the canonical variables of the 


Hamiltonian Eq. 81 


The condition of being tangent to ^ = 1 

implies cos[r2] = —/r, which we can write as cos[H] = 
(—1)'^, with M an integer. By using this we find that 
Eq. is just a reformulation of Eq. which is the 
equation defining ,dy of a fix-line. This equation ad¬ 
mits always a solution, in fact the condition = 1, which 
reads 


|A| A, 
A |A, 


(-ir = 1, 


(89) 


is just the condition of existence of a fix-line Eq. 66 
Therefore the point dx where 'S'' = S^(dx) is tangent 
to the horizontal line S' = 1 identifies a fix-line. As we 
have already seen this fix-line is parameterized by r, and 
the relation between r and the fix-line dx , dy (see Eq. 65) 
reads 


= (l-r) 


= 4r(l — r). 


(90) 


Eor these two values of dx,dy we find that '^'x = 'Py = 


0 (Eqs. 591, and the angles (px,'Py have to satisfy the 
condition (fx + = Mtt. The relation of ^ to r is 

^ = 3r2 - 4 t -h 1. 


(91) 


B. Stability of motion 


The stability of a trajectory is characterized by the 
level lines as follows: 

A particle motion defined by the initial condition 
((d£c)oj^%) is stable if—1 < S^'o < 1 and its level line 
{dx-, S''{dx)) crosses both the lines SP = \ and W = —1 or 
twice the same line SP = \ or SP = —1. Only in this case 
the motion is bounded by the periodicity of the coordinate 
VL. 

We next discuss the properties of stability as a function 

ofe 


1. Region 1 < ^ 


Observing part a) of Fig. 17 the red line / = C separates 


two classes of level curves, the curves below the red line 
go to —oo for dx approaching Appendix B discusses 
the geometrical properties of these level lines. It is clear 
that if 1/v^ < 1, i.e. for ^ > 1, the red curve that 
separates the two classes of level lines is located between 
S/ =\ and = 0 so that no particle can be stable. Any 
pair (da;, SP) chosen in ^ < da; and —\<SA <1 defines a 
level line that goes to da; —oo, 0 without crossing 
SA = ± 1 . 


g. Region 0 < ^ < 1 


This case is also illustrated in part a) of Fig. 17 The 


in pink in which every initial {{dx)o,S''o) evolves on a 
level line that is bounded. Specifically as the particle 
trajectory {dx,S'') cannot exceed —1<S'' < 1, it follows 
that the particle bounces back and forth on the lines 
S'' = — 1, and SP = \ keeping the invariant constant. 
Particles on the right of the level line tangent SP = \ 
cannot be stable. It appears evident that the maximum 
region of stability, marked by a region of thin red lines 
is determined by the level line tangent Xo SP =1. As 
discussed in the previous section, the point tangent to 
the line SP = 1 is a fix-line. This point is indicated as 
(d*)+ in part a) of Fig. 17 On the same picture we 


also find the point (da;)i, which marks the region of total 
stability. This means that for ^ < dx < {dx)i particles 
are stable for any initial phase H. For (d^,)! < dx < 
(d*)+ only particles which have S'' in the proper range 
(namely ft in the proper interval, see Fig. 22) are stable. 


In the laboratory frame, this means that only particles 
with (px + 2ipy in the same range can be stable. The point 
(d*)+ becomes the more “external” point of stability. We 
note that for ^ = 1 the point 'f/y/f, becomes I, and the 
pink dashed area goes to zero. 


3. Region —1/4 < ^ < 0 

The properties of the level lines for ^ < 0 are dis¬ 
cussed and summarized in Fig. of the Appendix B. 
A version of the same picture with the relevant quan¬ 
tities is shown in part b) of Fig. |17[ In this region it 
holds l/(2-v/^) > 1, that means that the blue curve 


level line tangent to = 1 determines a region dashed 


shown in part b) of the Fig 17 will always cross the line 
S'^ — 1 twice. The level lines on the left of (or above 
of) the blue curve exhibit a minimum and a maximum, 
while the level lines on the right of (or below of) the blue 
curve have only one maximum. We find that there are 
two curves tangent Xo S'' = 1: one above the blue curve 
tangent in (d*)_, and one below the blue curve tangent 
in (d*)+. In the dashed blue area we find the level lines 
with a minimum and a maximum: we call for convenience 
these lines of the “second order”. The minimum of these 
lines has coordinate dx always smaller than {dx)t, and 
the dynamics of particles in this region is to be bounded 
by the same curve S'' = 1. This means that for parti¬ 
cles in the dashed blue region their trajectory is always 
bounded in 0 < <1. The red dashed area charac¬ 

terizes instead the initial conditions where the level lines 
have one maximum, we call for convenience these lines 
of the “first order”. These first order lines always cross 
both lines SP =\, and S'^ = — 1, hence this region is sta¬ 
ble and as discussed in the previous section the level line 
tangent to the line SP = \ identifies a fix-line. Therefore 
(d*)+ is on the fix-line and characterizes the outer point 
of stability. 

The point {dx)i marks the region of total stability: if 
0 < da; < (da;)i particles are stable for any initial phase 

n. 












a) 

& 
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FIG. 17: Impact of the level lines on the stability of motion. 
The three pictures refer to the three relevant regimes. Note 
that in a) the curves stop at ^ because of the condition 
of Eq. 1^ 


In this regime, referring to Fig. [29lwe find that the 
maximum of the blue line is l/(2-\/—< 1, and the 
inflection point of the red line is l/(-\/—3^) >1. As a 
consequence all level lines below the blue curve describe 
unstable dynamics. Among the level lines in the region 
between the red and blue line, there are two level lines 
tangent to the line = \. The situation is illustrated in 
part c) of Fig. 17 One level line is tangent to = 1 in 
(o*)+, this curve also intercepts the line ^ = 1 in {ax)i- 
The second level line is tangent to = 1 in (a*)_. Only 
the blue dashed area is stable as any initial condition in 
the system of coordinates (dx, will move on a level line 
that crosses twice the line ^ = 1. Interestingly, the point 

V = ((a*)_, 1) is a point that does not move because it 
also identifies a fix-line (the level line is tangent to = 1 
in (d*)_), but at the same time the level line tangent to 

V has ^ > 1. Therefore any point in the neighborhood 

of V with — 1 < <1 will always remain close to it, 

and V is an equivalent of a fix point in this system of 
coordinates. We call this a “stationary” point. 

Lastly, we also find the point (dx) 2 , which is the ex¬ 
treme point of stability, for d^ < (dx )2 there is no stabil¬ 
ity. 


5. C < -1/3 

In this region the condition ^ < —1/3 implies 

1/-^—3^ < 1. This means that the inflection point of 
the red curve in Fig. is below or equal to line ^ = 1. 
Therefore if ^ < —1/3 no point in Oa, > 0, — 1 < < 1 

can be stable. 


C. Simulation examples 


In order to confirm the theory as developed in this 
paper we directly integrated the equation of motion of 
the coordinates and verified that the stability 

of particle is determined as previously discussed. The 
vicinity of the resonance is explored here on parameters 
that will be used later in multi-particle simulations. We 
use the following parameters A = —7.44 x 10“^ m“^, and 
Ar = 4.64 X 10“^, from which we find the scaling factor of 
Eq. 78 equal to A^/(16A^) = 0.024 m^, and the /i = —1. 


In Fig. we show an example of level lines for ^ = 
0.32, where each curve corresponds to a different /. The 
green curve is the level curve of the fix-line defined by 
T = 0.2, the corresponding invariant in the laboratory 
frame is C = 0.0155 m^. The set of orange dots in the 
picture shows the set of initial conditions that are stable. 
It is evident that the green curve bounds the set of stable 
initial conditions as discussed in the previous section. 

Part a) of Fig. 19 shows simulations in the regime 
— 1/4 < ^ < 0. The parameters of the plot are r = 











































a) 



FIG. 18: Level lines of the invariant. The green curve touches 
both the lines = ±1 and bound the stability region marked 
by the orange dots, which are the set of stable initial condition 
obtained from simulations. 



b) 



0.4,^ = —0.12. Among the level lines, two are of rele¬ 
vance and colored with green and red. The green curve 
is a level curve of the first order. This curve touches 
cos[n] = 1 in one point (fix-line), and the set of orange 
dots on the left of the green curve in —1 < ^ <1 are all 
the stable initial conditions. In the picture there is also 
a red line, which is of the second order, and this curve 
also touches cos[f2] = 1 in one point. 


19 the parameters of the level line 


In part b) of Fig. 
associated to the fix-line are r = 0.81,^ = —0.268, there¬ 
fore the level lines have pattern typical of the region 
— 1/3 < ^ < —1/4. The green curve is a level line of 
the second order and it bounds the stable initial condi¬ 
tions (orange dots) with the line cos [12 ] = 1. For sake of 
clarity we show in the part c) of Fig. 19 the same sce¬ 
nario in the da,, 12 plane. The role of the fix-line is very 
evident. It also becomes clear that there is a “stationary” 
fix point where the red level line touches the line 12 = 0, 
and a single red dot is in the center of an “island”. 



FIG. 19: a) r = 0.4, ^ = —0.12. Level curves of the invariant. 
The two curves correspondent to the fix-lines are in red and 
green. Note that the level line of the unstable fix-line (green) 
is bounding the stable particles; b) r = 0.81, ^ = —0.268. 
Level curves correspondent to the fix-lines, they are curves of 
the second order (see text). The orange dots are in the region 
of stable trajectories; c) the same situation as in part b) but 
now in the a^, 12 coordinates. 









































D. Summary 


Because of the resonance the motion of a particle is 
also affected by the time dependence of the variables 
ax,ay,ifxT'fy At any time the variables d^, dy are forced 
to lie on the straight line 2dx — dy = 2^, hence the motion 
of a particle is bounded, if dx remains bounded. The sta¬ 
bility properties are retrieved as follows: given the initial 
condition dx^dy^t^xi^fy, we find ^ and fl, from which 
is computed. The pair {dxT'3^) with ^ allows to identify 
to which region the motion belongs to and correspond¬ 


ingly to which type of stability (discussed in Sec. VIIIB 
Figs. fT7|). 


In order to classify the stability properties we will make 
use of two definitions of stability based on the discussion 
made in Sec. |VIIIB[ We call stability of type 1 the 
stability of a particle, when its motion is associated to a 
level line of the first order. Similarly we define a stabil¬ 
ity of type 2 when the stability is associated to a level 
line of the second order. 



The findings on the stability in the proximity of the 
third order coupled resonance Qx + 2Qy = N are sum¬ 
marized in Fig. 1^ In these pictures the fix-lines define 
the border of the domain of stability (pink line unstable 
fix-lines, and blue line stable fix-lines). As each fix-line is 
characterized by 0 < r < 1 , r parameterizes the border 
of stability in dx,dy coordinates as well. The relation of 
T with f is shown in Fig. 



T 


FIG. 20: (, as function of r. The curve shows a symmetry at 
T = 2/3, because for r > 1/3 the line 2dx — dy = C crosses 
the fix-line parabolic-like curve twice. 



FIG. 21: Summary of all the parameters that characterize 
the stability domain. On the bottom is a zoom for the sake 
of clarity. The border dividing the red dashed area from the 
yellow dashed area is the collection of the (oi,)! defined in 
Fig. 18 part a) and b). The green line is the collection of 
the (dx) 2 , i.e. is the part of initial condition that are stable 
beyond the fix-line. The blue line is the collection of the stable 
fix-line (a),)-, the pink line is the collection of the unstable 
fix-lines {&%)+. The points {dx)t are defined in part b) of 
Fig.fTTl 
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The stability properties are summarized in 
Fig. in terms of curves collecting the points 
{ax) 2 , (dj)-, (dx)+, {ax)i, {ax)t defined in Fig. The 
bottom picture of Fig. provides a zoom with the 
summary of all the symbols. We can therefore make the 
following general classihcation of the regions of stability: 

Complete stability The red wide area in Fig. is the 
region of complete stability, i.e. particles with axtiy 
inside this area are stable for any (px,<fy The upper 
edge of the red area is the set of the points (dj,)! 


identified in part a),b) of Fig. 17 


Partial stability type 1 The yellow dashed area is a 
region of partial stability that extend in the region 
— 1/4 < < 1. The allowed angles fl are given by 


\n - 7r(l + fi)/2 + 27rW| < AH 


(92) 


with N' an integer. The angle AH is represented 
in Fig. for the case of ^ > 0, and is obtained 


directly from Eq. 84 as 


Afl{ax,ay) = arcos 


1 Ifi - 


(93) 


with the invariant correspondent to the fix-line 
defined by ^ (or r). We here consider 0 < AH < tt. 



FIG. 22: Definition of AH for the case ^ > 0. A similar 
definition applies for other value of ^ seen in the different 
parts of Fig. 


Partial stability type 2 This region is found for 
— 1/3 <^<—1/4, and is shown by the blue dashed 
area in Fig. the stability is type 2, and H is 
bounded to oscillate around H = 0 or H = 27r. The 
green line is the collection of the points (ax )2 iden¬ 
tified in part c) of Fig. 17 and is the edge of the 
stability domain. This curve is not a collection of 


fix-lines, which indeed correspond to the blue curve 
(the stationary points). The stability angles AH 
are given by Eq. [Mj 


Partial stability type 1+2 In Fig. this region is 
the overlapping of the region of partial stability 
type 2 (blue dashed area) with the region of par¬ 
tial stability of type 1 (yellow dashed area). This 
area is characterized by a stability shown in part 
b) of Fig. f^for (aa:)i < cLx < {ax)t- In Fig. M 


the bend curve {ax)t is the collection of the points 
{dx)t defined in part b) of Fig. 17 In this region of 
Fig.[^ according to the initial condition a particle 
may have an evolution of H that is bounded to 0 or 
27r (stability type 2, blue dashed area of part b) of 
Fig.[^, or an unbounded evolution of H that spans 
periodically all [0, 27r] values (stability type 1, red 
dashed area of part b) of Fig. (Tt]). Also in this case 


the angles that are stable are given by Eq. 93 


With this theoretical considerations in mind we com¬ 
pare the stability domain as predicted by our theory with 


multi-particle simulations. In the part a) of the Fig. 23 


we show the stability domain as obtained from theory, 
in an dx,dy chart. The color code provides the informa¬ 
tion on the allowed angle AH. The picture retrieves the 
discussed properties and shows that the allowed range of 
angles becomes smaller when approaching the unstable 
fix-lines. The picture also shows the partial stable region 
beyond the stable fix-lines (blue line from theory). In the 
red area the motion is stable (complete stability). 

Part b) of the Fig. 23 shows the same result obtained 
by tracking. In order to test the theory we considered for 
this simulation a constant focusing lattice where a single 
harmonics was created. The strength of the driving term 
is A = —7.44 X 10“^ m“^, and A^ = 4.64 x 10“^. The 
stability domain was sampled with 100x100x51x51 initial 
conditions, and the stability was assessed over 1000 turns. 
The pink and blue curves in the picture are the collection 
of fix-lines as from the theory. These results are in good 
agreement with the theoretical predictions. 

Part c) of the Fig. 23 shows the stability domain ob¬ 
tained when a single sextupolar kick is used on a constant 
focusing lattice. In this case we are not in the ideal condi¬ 
tion of the theory, stated in Sect. IV, and some deviation 
is to be expected, resulting in an enlarging of the stabil¬ 
ity domain closes to the ax axis. For this simulation the 
parameters are A^ = 4.64 x 10“^, and K 2 = 2.52 x 10“® 
m“^ that makes A = —2.36 x 10“® m“^. 


Note that in spite of the fact that the strength of the 
sextupole is very small, part c) of Fig. 23 respects the 
scaling expected from the theory (apart from the devia¬ 
tion along dx)- 


For more realistic cases in which a lattice with alter¬ 
nating magnet is used, or an arbitrary distribution of 
sextupolar errors is employed, this scaling will be similar. 
This subject will be discussed in the following chapters. 
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FIG. 23: Complete stability domain for a constant focusing 
circular accelerator. Part a) the stability domain as obtained 
by the analytic theory; Part b) the same result obtained by 
a particle tracking, with a computer code, of 100x100x51x51 
initial conditions, and particles are tracked for 1000 turns. 
Part c) shows the stability domain obtained for one sextupolar 
kick (to be explained in the following chapter). The result 
shows an enlarging of the stability domain closes to the dx 


IX. ANALYSIS IN AN AG CIRCULAR 
ACCELERATOR 

Previously we have made the analysis for the constant 
focusing case for a resonance excited by a single sextupole 
term, we will repeat here the analysis in the case of an al¬ 
ternating gradient (AG) structure. We discuss the effect 
of a distribution of errors and use the strategy previously 
discussed to formulate a theory of particle motion and 
fix-lines stability. We study the onset of stable motion in 
the proximity of a third order coupled resonance. 


A. AG structure set-up 


We consider the equations of motion in an AG struc¬ 
ture. The equations of motion take the form 


dPx 

ds"^ 

dh 

ds"^ 


+ K{s)x = fx{s,x,y), 
+ kv{s)y = fy{s,x,y). 


( 94 ) 


Here the kx{s),ky{s) are the “time dependent” focus- 
ing/defocusing quadrupole strength. The solutions of 
these equations for fx = fy = 0 are 


x{s) = \/l3xax cos{<px{s) + (fix), 

y{s) = y/Pyay COS{(t)y{s) -h <fy), 


( 95 ) 


where /3x,l3y are the well known beta function and 
'/'a:(s) = Jo/3x(s)~^ds, 4>y(s) = j3y[s)~^ds. The 

Hamiltonian of the system has the form 


Hn = 


{x'? , , , iy'? 

— + kx{s)- + — 


+ ky{s) — , ( 96 ) 


where x' = dx/ds,y' = dy/ds, and the canonical equa¬ 
tions are 


dx 

dHo 

dx' 

dHo 

ds 

dx' ’ 

ds 

dx 

dy 

dHo 

dy' 

dHo 

ds 

dy' ’ 

ds 

dy 


( 97 ) 


B. Equations of the constants 


Despite the fact that the unperturbed motion is time 
dependent, the equations of the otherwise constants 
ttxT^Px-iO.yi 'fy for the perturbed system with Hamiltonian 
H = Hq + Hi are as before 


dax 


dipx 

=2 

ds 

difx 

ds 


day 


dify 

=2 

ds 

dify 

ds 



dHi 
dux ’ 
dHi 


axis. 


( 98 ) 
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where we set the perturbing Hamiltonian Hi due to dis¬ 
tributed sextupoles as 


Hi 




(99) 


Let’s consider a distribution of sextupolar errors localized 
at specific positions Sj. The function ^ 2 ( 3 ) in this case 
reads 


To this purpose we call 

g 

g 

and Pxj = Pxisj),l3yj = Pyisj). Hence 


(106) 


k 2 {s) ='^2K2jS{s — Sj + mL), (100) 


where L is the periodicity of the system, hence the length 
of the accelerator, and K 2 j is the integrated strength 
of the sextupolar error located at the position Sj. The 
unperturbed solution reads 


X = \/l3xax 


_^iAx _|_ g iAx 


y = '/PyO-v 


_ ^iAy _j_ g 'i'Ay 


( 101 ) 


with Ax = 4>x{s) + ifx, Ay = 4>y[s) (fiy. We substitute 
(x, y) into the Hamiltonian Hi and obtain 

Hi ='^K2jS{s - Sj + 

-^K2jS{s - Sj + mL)^^yPxaxl3yayX 


(e ^ +e ")(e ’' + e “) . 


( 102 ) 


Note that the function (j)x{s) has the property 

(j)x{mL) = 2TrQxm, (103) 


Hi = ^K 2 jS{s - Sj +mL)^{l3xjaxf^'^ie^^^ +e 

3,HI. 

- ^2 H 2 jS{s - Sj + mL)^^JPxjaxPyj x 


3,HI 

T-y{ 


fly(e*^"^-I-e -i^e 


(107) 


We substitute the delta expansion and Axj,Ayj obtain¬ 
ing 




3 / 2 / 3 


j t 9=0 


exp{i27r^—-I- i{&x{sj) + 2 ttQxy + ‘fix){2q — 3]} — 

1j 1j 


161/ \/ Py:j Pv 


j t q—o p=0 


1 \ / 2 


g \p 


X exp{i[27r—-|- &x{sj){2q — 1) -f %{sj){2p — 2)]}x 

X exp{i27r[Qi:(2(} — 1) -I- Qy{2p — 2) -|- t]y}x 

Jj 

X exp {i[ipx{2q - 1) -b Pyi^p - 2)]}. 

(108) 


with m an integer {(px is defined for the linear case). 
Therefore we can define the function 

^x{s) = (t>x{s) - 2 ttQx^, (104) 

which is a periodic function with ^xi'm-L) =0. In a 
similar way we define ^y{s). Therefore 

g 

Ax — ^x(^) ~b 2.TTQx~r “b y^x^ 

S (105) 

Ay - ^y{s) 2TrQy— ify. 


We therefore find that the Hamiltonian Eq. |102| is com¬ 
posed of the periodic functions ^^(s), ^y{s), I3x{s), I3y{s), 
which have the periodicity of the accelerator L, and 
also by non-periodic functions like the term 2'!TQxsjL 


(etc.). On the other hand the delta function in Eq. 102 


is different from zero only at the longitudinal posi¬ 
tions Sj — niL, and at these positions the periodic 
functions ^x{s), S’y{s), f3x{s), Py{s) assume the value 
^^x(sj), {sj), f3x{sj), f3y (sj) independently on m. We 
can therefore replace in all the periodic functions of 
Eq. 102 the coordinate s by sj. 


As the tunes Qx,Qy are close the resonance, meaning 
Ar = Qx + 2(5y — is small , it follows that between all 
oscillating terms in Eq. |108| only two have a very slow 
frequency of oscillation, while all the other can be con¬ 
sidered as fast oscillating. We neglect the fast oscillating 
terms and keep just the slowly oscillating ones, that is we 
select g,p, t in Qx{2q—l) + Qy{2p—2)+t such as to obtain 
slow frequencies. This happens for g = l,p = 2, t = —N, 
and q = 0,p — OA = H . Consequently the slowly varying 
Hamiltonian reads 


~ op '/oAo.y 'y ( K2j ijPxjPy 


COs[27r^A -I- ^x(sj) + 2&y{sj) -I- 27rAry + ipx+ 2q}y\, 
Lj 1^ 

(109) 


we re-write Hgi in the following form 


Hsi = i/a^CLy { Ac cos (fix + 2vJj, -b 27r A. 


Li 


—Ag sin 


y^x “b 2.py -\- 27 tA x } 

1j j 


(110) 
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with 


Ac — ^ Ar2j\//3a;jA 


iX 


X cos j^ 27 r^iV + ^xisj) + 2 ^y{sj) 


As — \/Pxj Pyj X 


( 111 ) 


3 

X sin 


s,- 


2^^7V + ^,(s,) + 2%(s,) 


The coefficients Ac, Ag depend only on the distribution of 
the sextupolar errors, on the optics, and on the selected 
resonance (in this case Qx + 2Qy = N). Next we remove 
the time dependence in the argument of the trigonometric 


function in Eqs. 110 with the transformation 


-^x — ^X 


Qjy - Qjy 


^x — tx^^iTT^ 


( 112 ) 


ty^TT 


where we require that 


tx + 2ty — 1. 


(113) 


Hsi =y/axay {Ac cos [(px + 2(py] - Ag sin [(px + 2(^y]} + 

/- ~ ^r, Ac 

~\~\^Ctx^x 4“ )27r . 

(114) 

The last term is equivalent to the similar term in the 
Hamiltonian for the constant focusing Eq. It arises 
from the requirement of making the Hamiltonian time 
independent. At this point it is convenient to define the 
angle a as 


The new Hamiltonian reads 


Ac 


Ag 


cos a = sm a = — 

A A 


(115) 


with A = a^A^ + A^ > 0. Therefore the final form of the 
slowly varying Hamiltonian in an AG lattice becomes 

27rAc 
(116) 


Hgl — A-y UxCiy COs(^(px 2 ipy cr) -t- {jXx^x 4 “ 


The canonical equations are 

-a'x = 2 = 2^/h^dyAs\n{^px 4- 2(py + a) 

Cf'^x 

-! c 1 ~ A f, - A 27rAc 

tPx = 2 = 2 ay Acos{(px + 2ipy + a) + G— 


dda 

dHsi 


= A^/^dyAsin{(px + 2(py + a) 


®y ~ ^ pi~ 
d(py 

~f rt /-—A /- ~ \ 27rAc 

P>y = = 2y'axAcos(ipx + 2ipy + a) + ty ^ 


(117) 


multiplying the first equation by 2 and subtracting the 
third one we find that 


2(3!x ^y — 


(118) 


where C is a constant. The previous expression is then 
an invariant of motion. 


C. Special solutions 


Via Eq. |117| the equations of the fix-line are 

—d'x = 0 = A2\/^dy sin((,5a; -|- 2(py + a) 

, A 1 / 27rAc 

(Px = 0 = A2 ^ p^ fiy cos[(px + 2ipy -\- a) +tx 


. ^ - y - -^y-rx I -T-y ' —/ ' ~X y 

V ^x ^ 

-dy = 0 = Adi/^Uy sin((^c -|- 2(py + a) 

<p'y = 0 = A2^/d~x cos{(px + 2(py + a) + ty 


The condition a(, = = 0 is obtained from the condition 


which yields 


sin(i^2: 4- 2ipy -I- a) = 0 


ipx + 2<py + a = ttM 


(119) 
dition 

( 120 ) 


( 121 ) 


with M an integer. That means that the other two equa¬ 
tions become 


0 = A2^dyi-l)^ +tx^^ 


2'\/ Ua; 

0 = A2v^(-l)'^+t, 


271 Ar 


( 122 ) 


Multiplying the second by 2 and summing them up and 
recalling Eq. |113| we find that the equation for the fix-line 
reads 


0 = A 


1 


-a. 


y 


(-1)" 4- 


27rAc 


(123) 


/ax " / L 

and the condition for its existence is that 

AAc(-l)“ < 0. (124) 

Therefore from the equations for the fix-line we write 


0 = A- 


1 


(-l)"Ac + 


27rA^ 


L 




0 = A2\/d;(-l)'^Ac-f 


(125) 


from which we conclude that tx > 0 and ty > 0. That 
means that G > 0 and 1 — tx > 0, so that we get 1 > 
tx > 0. We then parameterize tx,ty as follows 


tx 

ty = x(l - t) 


(126) 
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with 0 < T < 1. We write the solution for the fix-line as 
a function of r, and get 


ay — 


(27rAr)^ 

4A2L2 


t(1 


r), 




(27rAr)^ 

16A2L2 


(l-rf. 


(127) 


We find the same solution as for the case of the constant 
focusing except that the scaling depends on A/L, with 
A being the driving term of the resonance as already 
computed by Hagedorn/Schoch [SHS], Guignard [6l[7]. 


D. Fix-line 

The analytic form of the fix-line can now be found. 
From (p'^ = (p'y = Q we get 


and returning back to the coordinates of the Hamiltonian 
Hsi we find 


g 

— ^y,o ty27rAj. —, 


(129) 


with sin((^a;_o + 2(^y_o + c«) = 0- The solution has therefore 
the form 


X = \/fixax cos |^0 x(s) 

tx‘^'^Ar “h (^x,oj 

/ ,— r j 

x('5) tx27vAr~ ^a;,oj ’ 

- \/^xax sin [(/>x(s) 


y = \/?yay cos ^(l>y{s) 

iy27rAj"— 

y = /^V^cos U 

V Py 

tiy27rA)"— 


i Stt 



— Pyay sin — ty27rAr|- + CPy,oj - ty27vAr^j 

(130) 

This is the equation of the fix-line, which is a closed line 
in the 4D phase space. 

The Poincare surface of section is identified by the con¬ 
dition s = LAf with Af an integer correspondent to the 
Af-th turn. From the solution of the fix-line equation we 
have sin((^a;_o + + a) = 0 and we find 


(px,o + -I- a = ttM, (131) 


with M an integer properly taken to guarantee the exis¬ 
tence of the solution. It is straightforward to show that 

{Qx ~\~ ^x,oA~ 

2(Qy - tyAr)2TrAJ' + 2ify^Q = (132) 

2'KNJ\f ~\- ttM — a. 


Therefore, by using Eq. m we hnd that on the 
Poincare surface of section we have 

X = \/Pxttx cos [—2{Qy — tyAr)2'Kj\f — 2ipy^o — a + ttM] , 
V — \/^y^y [(^y tyAjA^2TTj\f 7^y,o] •> 

(133) 

with (py^Q an initial phase that sets the starting point. 
Lastly, we find the result that the fix-line can be param¬ 
eterized as 


Xt — A^x^x cos[ 2{C^y tyAx^t (X -f 7rAf], 
yt — '\/l^y^y COs[(t^y tyA,.)t]. 


(134) 


The presence of the angle a is the only difference with 
respect to the case of the constant focusing lattice. 
However, by changing the longitudinal position of the 
Poincare surface of section, a can be made equal to zero. 


E. Stability 

The discussion of the stability of the fix-line is the same 
as for the case of the constant focusing system. The 
conclusion holds as before: for r > 2/3 the fix-line is 
stable, and for 2/3 < r < 1 the fix-line is unstable. 


F. Stability of Motion 

Next we study the onset of unstable motion. From the 
canonical equation Eqs. m we keep the first equation; 
due to 2ax — Qy = C the third equation is redundant. 
The second and fourth equations can be summed up to 
yield 


—a'x = 2-\/^dyAsin(r2 -|- a) 


n' = 


‘^'s/ ^x 


with Q = (px + 2(py. 
equations as 



Acos(r2 -I- a) 


27rAf. 


L 

(135) 


Therefore we re-write the set of 


—a'x; = 2^/^{2dx — C')Asin(H -|- a) 


n' = 


—T={2dx — C) + 

V 


Acos(H -I- Of) 4 


The function 


I{dx,^) = 2Ky/^{2dx — C) cos(H4q;)40; 


27rA^ 


27rAj. 

L 

(136) 

(137) 


is the invariant of the system as obtained for the constant 
focusing discussion, except for the presence of the angle 
a in the cosine function, but this is just a shift of the 
stability diagram. Hence, the conclusions obtained for 
the case of the constant focusing structure hold in the 
AG structure as well. 

























26 


G. Comparison with simulations 


In Fig. we compare the reconstruction of the stabil¬ 
ity domain with part a) of Fig. 25 that is derived from the 
theory and numerical simulations. As a practical exam¬ 
ple part b) of Fig. 25 shows instead the same result for the 
SIS 18 reference lattice equipped with a single sextupolar 
kick. As for the case of the constant focusing we find 
an enlarging of the stability near the axis. Part c) of 
Fig.j^shows the stability domain for the SIS18 where 12 
sextupolar errors are randomly excited. The sextupoles 
are all placed at the beginning of each super-period. In¬ 
terestingly the results are still close to those predicted 
by the theory. A small shrinking of the stability is found 
anywhere close to the axis Ux- The discrepancy can be 
attributed to the effect of high order harmonics not in¬ 
cluded in the theory. 
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FIG. 24: Comparison of the wavelength from analytic theory 
(blue dots) and retrieved from simulations (horizontal color 
black and vertical color red). This result is obtained for a 
FODO cell with one sextupolar kick. The agreement is ex¬ 
cellent with a discrepancy that increases in the proximity of 
T = 2/3. 


We have verified this ansatz by simulations in a FODO 
structure where we artificially made a single harmonics 
stronger by placing a sequence of sextupoles with alter¬ 
nating sign of strength. The results are shown in part d) 
of Fig. |25| the stability domain very well retrieves the 


theoretical one of part a) of Fig. 25 


Lastly, we tested the prediction of the secondary tunes 
for a particle slightly off a fix-line. These results are 
shown in Fig. The simulation is made in a FODO 
cell equipped with a single sextupole. The secondary fre¬ 
quency is measured by comparing the side bands to the 
nominal tunes. The theoretical prediction is obtained 
from Eq. where is replaced by 27rAr/L (green 
dots). From the simulations we retrieved the wavelength 
in x,y (red dots). The picture shows an excellent agree¬ 
ment except for slight deviations of the side bands where 
the curves diverges and high order terms become relevant 
(near r ~ 2/3). 


a) An b) An 
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FIG. 25: Gomplete stability domain for an AG circular ac¬ 
celerator. Part a) the stability domain as obtained by the 
analytic theory for the SIS18 lattice with a distribution of sex¬ 
tupolar errors. Part b) shows the stability domain of SIS18 
obtained for one sextupolar kick. The result shows a slight 
enlarging of the stability domain close to the axis. Part c) 
the same result obtained by tracking of 100x100x51x51 initial 
conditions, and tracked for 1000 turns. Part d) is an example 
of a FODO cell with a distribution of sextupoles such as to 
excite only one harmonic. We find a good confirmation of the 
theory at the edge of the stability domain. 


X. TORI CUT 


In this section we discuss the Poincare surface of sec¬ 
tion of the tori surrounding the fix-line. We have already 
obtained those orbits from the numerical simulations in 
Figs. [T^ We show here that these orbits will be 

retrieved from the invariants. In those figures the orbits 
in X — x' are obtained by tracking one single particle for 
many turns, and plotting x,x' at each turn when the 
condition y'/y < e and y > 0 are satisfied; here e is a 
positive small number. The orbits in y — y' are obtained 
by making a selection of the particle when x'/x < e and 
a: > 0. 

Without losing generality we make the discussion for 
the case of the constant focusing structure as presented 
in Sect. IV, V, VI, VII, VIII. Our starting point is to 
consider a tori, identified by the invariants I = I{ax, O), 
and 2ax — o,y = C, filled with particles. That means 
that a particle in the tori have ax,ay,^ satisfying the 
two invariants, with phases (px, ‘Py satisfying tpx -f 2(py = 
O (the fl variable is defined in Eq. 70). The situation 
is equivalent to what is done in the simulations when 
tracking one particle repeatedly for many turns, in the 
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condition that the particle explores all the tori. 

The spatial coordinates of the particles populating the 
tori are given by 


‘normailzed” coordinates 


— COs('0a;), 


y 


,cos{tjjy) (138) 


VFx " 

where we dehne for convenience 

'^x — Qx^ y^X} '^y — Qy^ “t” 

where the quantities ax,ay,ipx, Vy vary in time. For this 
case (constant focusing) the beta functions are j3x = 
1/Qx,l3y = l/Qy As discussed before, a specific set of 
initial particle coordinates is identified by the initial con¬ 
dition axo,ayo, ‘PxO,‘PyO! and this condition is used for in¬ 
tegrating the equation of motions of (these 

are the Eqs. 35 with J^x = -Fy = 0). The momentum 


coordinates are obtained by differentiation 



^\J ^x 


cos(V’x) - '/a^SlTL{'ll)x){Qx + ^p'x), 


cos('0y) - Ja:ySu\{'ij}y){Qy + ip' ). 


(140) 

The variables a'x,a'y,p'x, p'y can be expressed in terms of 
a'xTay,p'x,p'y, by using the transformations 


Px - Px ^x 


dy — 

= ^y 


— ^1! ^y^ ■ 


(141) 


and in these variables we know what is a'x,a'p'x, pL from 


'^y') 'r a: 5 ry 

the equation 35 (with Jq, = = 0). By substituting the 

equations of motion Eqs. [35] in Eqs. |140[ and after some 
algebraic calculations we obtain 



= kcLy s\tl[-iI)x + VL) - yjhx shi{il)x)Q. 


= 2KyJaxay sin(—■i/jy -|- 11) — y^sin(i/)y)(3y. 

(142) 

At this point we define the “normalized” coordinates 


(suggested by Eq. 78) 


^ = 


Kf 


y = 


Vp = 


X 

4 A 


Ar 

(d 

X 

ydO 

y 

4 A 

\/%j 

Ar 


y 


4 A 


= 


4 A 


(143) 


de^y 


4 A 


A. 


= y'\/% 


4 A 

A. 


Therefore the particle in the 
reads 

^ = \/%x cos{tjjx) 


Kp = pAr&y siTi{-'il>x + H) - \/^sm{'il)x) 

y = ^/hyCOS{ljPy) 


(145) 


yp = -y A^Ar v^a^^Oy sin(-'0y -|- H) - 


,sin(V>y) 


Oy are those dehned in Eq. 78 


where the coordinates 

Using the definition px + 2py = U, we find that ipxP^' 
are bounded by the relation ipx + 2ipy = NO -|- U. As the 
Tori is being hlled, we make sure which particles satisfy 
some special condition at say 0 = 0 (we do not need 
to consider many turns, because the tori is completely 
populated with particles, and after each turn the tori 
remain populated in the same way). Hence the tori cut 
at 0 = 0 yields 


'^X 2'i/?y - H. 


(146) 


In order to obtain the orbits in x — x', i.e. in ?c — ?Cp, we 
require that we select only the particles with yp = 0. For 
small Ar this happens for all particles in the Tori, which 
have ■i/’y = Trm, with m an integer. That means that we 
choose particles with a phase in the x plane as 


ipx = —27rm + H. 


(147) 


By substituting Eq. 147 into the Eq. 145 (formula for 
?(p) we find that the coordinates of the orbits obtained by 
cutting the tori imposing yp = 0 are given by the map 


K. = cos(H) 
!(p = -\/^sin(H) 


(148) 


where dx-,^ must satisfy the invariants that defines the 
tori. As we started with the filled tori, any pair da;,H 
satisfying / = I{dx-,^),2dx — dy = 2^ will identify a 
particle initially populating the tori, and can be used 
in Eq. 148 The condition y > 0, i.e. y > 0, implies 


m = 2n with n an integer. However, this condition does 
not change Eq. |148| 

We proceed in a similar way for the cut in = 0, we 
select the particles that satisfies ijjx = with m an 
integer, in addition the condition that cc > 0, i.e., a: > 0, 
implies m = 2n, with n an integer. As discussed before 
we consider the tori at 0 = 0. Therefore once we fix O, 
there are the following phases possible 


ipy = — — Trn. 


(149) 


Substituting it in Eq. 145 (last two equations) we find 
the following map 

n 


y = (-l)”VayCOS - 


yp = -(-l)”\/aysin ( ^ 


(144) 


(150) 
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FIG. 26: Tori cut for parameters Ar,A such that jj. = —1. 
The solid curves are the orbits obtained from the analytic the¬ 
ory, while the dots are the coordinates of one particle tracked 
and plot when it satisfies the tori cut condition. The red 
curves are the orbits of particles in the red dashed region of 
part b) of Fig. (stability type 1); The blue orbits are for 
particles with initial condition in the blue dashed region of 
part b) of Fig. |17| (stability type 2); Note in the center of the 
circular blue orbits one red dot: this is the “stationary point” 
discussed in Sec. IVIIIB4I 


-(- 1 )^ 


for fj, = 1, and a,t y = (— 


% 


= 0 for 


fi = —1 under the tori cut conditions. Figure 26 shows 
a comparison of the orbits in a; — a::' and y — x' with 
the particle tracking. For this simulation it holds y = 
— 1. The color of the curves is associated to the stability 
properties of the level lines as discussed in Fig. [T^ 


XI. CONCLUSION 


We have derived the theory of the fix-lines and the 
stability domain close to a third order coupled resonance. 

We find that the number of fix-lines is infinite, and 
that they are of stable or unstable type. The set of fix¬ 
lines form a continuous curve of parabolic shape in the 
plane of the single particle emittances d^, Oy. The theory 
suggests a natural parametrization of the fix-lines based 
on the canonical transformation that makes the system 
time independent. 

Our theory predicts the behavior of perturbation close 
to a fix-line and we find the property that the perturba¬ 
tion oscillates along a line with 2Sax = Sdy. Our analysis 
shows that these objects define the border of stability 
when a single harmonics is excited. The structure of the 
canonical equation of the constants suggests a natural 
symmetry set by the resonance condition. We find that 
for the variables the dynamics imposes an invari¬ 

ant, and using this invariant we discuss the complete set 
of topological properties of the level lines. It turns out 
that the level line associated to the fix-line creates a bar¬ 
rier to particles, hence they set the edge of the stability 
domain. 

The extension of the theory to the alternating gradi¬ 
ent circular structure has been worked out essentially re¬ 
trieving the driving term of Hagedorn/Schoch [SHS] Guig- 
nard 013 , setting then the relation between the fix-lines 
and the onset of unbounded motion. The distributions 
of errors creates mainly an angle a that enters into the 
theory, but it does not otherwise change the conclusion 
found for the constant focusing case. 

This work will be the foundation for the study of the 
resonance crossing induced by space charge conjoint with 
a coupled resonance line. 


As before dx, 12 must satisfy the invariants that define 
the tori, i.e. I = I{dx, 12), 2.dx — dy = 2^. We recall that 
the fix-lines are defined by the level lines of / = I{dx, 12) 
that are tangent to the line ^ = 1, with '3^ = —p, cos 12, 
being ^ = ±1, the sign depends on the combination of the 
signs of A, Ar- That simply means that 12 of the fix-line is 
given by 12 = 12j,/ = |(1-|-^). Substituting this value into 
Eq. 150 we find that the fix-line are located a,t y = 0,yp = 
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(151) 


XIII. APPENDIX A: THE PERTURBED 
MOTION AROUND A FIX-LINE 

The starting point are the equations of the perturbed 
system, Eq. that we report here for convenience 

5a'^ = —XSifix — ‘2XS(fy, 

Sd'y = —2\5(px — A\5(py^ 

-^xy^^y) 

^^y ~ -^yx^dx -^yy^^y^ 

after some algebra they become 

5d^ = UJ Sdx T X(^Axy 
6dy = — uj'^Sdy — X(Axx + ‘^Ayx)Cp, 
dipx — X(Aj;x 2Axy^{S(^x ~X 2(5(py), 

S<fy = — X{Ayx + 2Ayy)[6iPx T 2(5(^y), 

where we defined 

= X{Ax;x + ‘^Ayx + 2.Axy + '^Ayy). 


with integration constants Bx,By,Tx,Ty. The equa- 


(152) 


(153) 


The coefficient Cp is a constant obtained from integrating 
Eqs. |151| 


26dx = Sdy + Cp. (154) 

From equation Eq. |152| we also find 

{Sifx + ‘^5(py)" = -uj'^{5(fix + 2(5(^y). (155) 

The general conclusion is that for 

> 0 (156) 

the quantity S(fx + 2Sipy is stable and oscillate around 
zero as well as Sdx, and Sdy . We now consider the case 

we find 


uj > 0. Integrating Eq. 155 


S(fx + 26ify = A cos(a;0 + ^) 


with A,^ some constant. Using this result in Eq. |152| we 
find 


Sifx — — X{Axx + 2 ^ 3 ;^) A cos(w0 + ^), 
S^y — — X(Ayx 2Ayy) A cos(cu^ T ^). 

The integration of these equations yields 

A 


(158) 


S^Px — X{Axx + ^Aj;^)-^ cos(u;0 + ^) + CxO + Dx, 

S(Py =X(yAyx + 2Ayy) -^ COS^LUO T ^) T CyO T Dy, 


(159) 

with Cx,Cy, Dx, Dy the integration constants to be de¬ 


termined later. The first two equations of Eq 152 can be 
integrated as well, we find 

Sdx —BxAivijjjO -\- Tx') 2 ^Axy ‘^Ayy^Cp, 

(160) 

Sdy =ByAYi{i09 ~\- Ty'j ^(A^ja; “t" 2Ayx^Cp, 


tions |159| and 160 are obtained from Eq. 


solutions must be consistent with Eq. |151 
the constraint imposed by this condition. 

We use the first two equations of Eq. |151| 

Sd'x = —XS(px — 2XS(py, 

Sd'y = —2XS(px — dXSpy, 

and substitute Eq. |160| and Eq. |157[ We find 

BxUi cos{ui9 + Tx) = —XAcos{(jj9 -I- ^) 
ByUj cos(uj9 -I- Ty) = —2 XAcos{(jj9 + 


152 but these 
We next find 


(161) 


(162) 


therefore it must be = Ty = = —XA/uj^ and 

By = —2XAju. The solution becomes 

Sqix — sin(w^ T T 2 
Sciy = sin(u^^ T 2 

A 


—X{Axx + cos(u;^ + ?) T 

—^^-^yx + ^^yy) —2 T ^) + CyO + Dy^ 


Now we use the last two equations of Eqs. |151 

SiPx - Axx^^X T AxySOiy 

- AyxScix T Ayyddy,^ 

and substituting the Eqs. |163| we find 

Cx 2 ^Cpi^AxxAyy AxyAyx')-, 


(163) 

(164) 


UJ 


Cy o^Cpi^AyyAxx -^yx-^xy')’ 


(165) 


UJ- 


(157) Therefore the solutions take the form 


Sdx — -sin(a;^ + ^) H- -(^Axy + 2Ayy)Cp^ 

5dy = sin(u^^ T ^) i^i^Axx T ‘^Ayx^Cp^ 

UJ UJ‘^ 

^^x —X{^Axx “t“ ‘2Axy^ 2 cos(cu’^ “h ^) 

~\~ 2 '7^Cp(^AxxAyy AxyAyx^O T Dx-, 


(166) 


UJ 


A 


S(Py - X(^Ayx “h 2j4yy) ^ COS(^Uj0 T ^) 


^yx 

X 


2 Cp(^AxxAyy AxyAyx'jO T Dy. 


As 


5(fx + 2S(fy = Acos(a;^ + ^) 
substituting we find 

Dx + 2T)y = 0 


(167) 

(168) 
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Therefore the general evolution of a perturbation around 
a fix-line point takes the form 


— sin(uJ^ T T 2 ^-^xy T 
^Q'y — sin(u^^ T '^) o^-^xx T ‘^-^yx)^pi 

LJ LJ^ 

^^x — ^(^-^xx “t“ ‘^^xy') 2 COs(ctJ^ “h 

+ 2— ^Cp^Axx-^yy ~ AxyAyx)9 — 2Dy^ 


UJ- 


A 


5(py —X{Ayx + 2Ayy)— cos{uj9 H- 


A 


9 ^p{_AxxAyy AxyAyx^9 “h Dy. 


UJ 


(169) 


The free parameters A^ (7p, Dy are found from Eq. 169 
at ^ = 0 using the initial conditions 


{Sq.x)o — 2 

UJ UJ 

{^^y)o ~ T ^Ayx^Cp, 

UJ UJ 

{_9(Px)o —Xi^AxX “t“ 2y43;y) — COS(^) 2Dy^ 

- X(^Ayx 2Ayy'j — COS(^) “h Dy. 


(170) 


A. Discussion 

Suppose we consider a fix-line defined by d^, dy This 
means that a particle on this fix-line has constant d^^, dj,, 
while its phases (fx, ‘fy are also constant, and satisfy the 
condition ipx + 2(^y = ttM. Taking any other pair of 
ifx, <f>y that satisfy (px + ‘^'Py = ’’"Tf simply means to take 
another point on the fix-line, but again (px,Py remain 
constant. Therefore any phases (px,Py satisfying (px + 
2(py = ttM are equally good. 

By saying that ax, ay identify a fix-line we intend the 
identification of the fix-line without any reference to the 
phases of a particle on it. 


Consider now at 9 = 0 a point of coordinates 
dx,ay,(px,Py lying on the fix-line identified by dx,dy. 
Consider now a perturbation 5dx, Sdy, Spx, Spy from that 
point, and consider the particle with coordinates 


ax T ^ax, ay -t- Say, px Spx, Py Spy. (171) 


The evolution of the coordinates of the particle is given 
by dx = const, dy = const, (px = const, (py = const, and 
by the evolution of the perturbation Sdx,Sdy,5px,S(py, 
which is given by the equations Eq. |169[ We note that 
if Cp = 2(dda;)o ~ {^^y)o ^ 0) the oscillation of the per¬ 
turbation Sdx, Sdy is not around dx,dy, while the pertur¬ 
bations S(px,Spy are unbounded. The perturbative anal¬ 
ysis seems to col lapse . On the other hand the quantity 
Spx + 2Spy (Eq. 1671 remains bounded around zero! 

Lets now consider the point with coordinates 


Ox —Ox T SOx — Ox 2 ^-^^y 2^yy)Cp, 

dy —dy -f Sdy = dy '^^{S^xx T ‘2jA.yx)Cp, (172) 

pI =Px + ^pI =Px- 2a, 
pI =Py + ^pI = Py + ot, 

where we defined for convenience Sd%, Sdy, Sp^, Spy. The 
parameter a is an arbitrary parameter that defines the 
point. We observe that the coordinates dx + Sdx,dy + Sdy 
oscillate around , aj). 

By construction the coordinates, d%,dy,p%,Py do not 
change in time. We will prove that they identify a fix-line 
point of a new fix-line identified by d%,dy. 

We start by writing the equation of the fix-line in the 
general form 


dHsi 2^^ 

ddx ddy 


= 0 , 


(173) 


this equation is Eq. 42 (top) expressed in terms of the 
slowly varying Hamiltonian Hgi- We no w check if a°, dy 
satisfies the previous equation (Eq. |173 ). We find 


dH. 


si 


dH. 


(~C ~C\ I /~C " C \ '-'-^-*•31 / — ~ \ I -^-*-31 / — ~ \ f~c i 

[Ox, Oy) -f 2,—^ \Px^ 1 \^x, ay) H [ax, ay)6ax + 

y 


dHsi 


d^Hsi 


d'^Hsi 


ddx “ dd„ ddx ddi ddxddS^' y 


[ax, S.y)Sdy-\- 


dH,,i ~ , r, di^Hsi ~ 

+2 (ca;, Oj,) -|- 2^_ [ax,ay)5ax + 2[ax,ay)Say. 


(174) 


ddyddx 


ddyddy 


Taking into account that dx, dy is a fix-line, then 
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dHsl,~c ~c\ , ~c ^ d'^Hsi . 

■(®a:)®y) + 2 [a^,ay,)— {ax,ciy)oa^ + (a^;, 0 ^) 00 ^,+ 


dux 


+2 


dal 
d^Hsi 
ddyddx 


ddxddy 


0^ H 

(^Ux ■> 2 )^(2y. 


(175) 


aa2 


Using the definitions for Axx, Axy, Ayx, Ayy the previous 
equation becomes 


da si ,~Q ~c\ /~c ~C \ 

-^iax,ay) + 2-^iax,ay,) = 

+ 2Aj,2;)5a^ + {Axy + 2Ayy)5dy\. 


(176) 


Now we substitute 5d%, Sdy from Eq. 


172 


and we find 


perturbation with respect to the new point reads 
5dx = - sin(a;'6» + ^'), 

LJ 

= ~ ^ sin(a;'6> + 

Sifx =^'i^'xx + ‘^^'xy)^ cos(ui'0 + ^') — 2Dy, 
5(fy =\'[Aly + 2A'yy)—^ cos(w'^ + ^') ^ H^, 


(181) 


dH. 


si 


1 A 


f)J-T 

{di,dy) + 2^{d'xra:,) = 


2^^p[i^xx + ‘^-^yx){-Axy + ‘^^yy) 

~ i^xy + ‘2Ayy)(AxX + ‘^Ayx)] = 0. 


(177) 


We conclude that a^,a° satisfy the equation for the fix¬ 
line, and therefore it identifies a fix-line. The phases 
also satisfy the condition of a fix-line, in fact we 

find 


ip'x + 2(fy = if>x-2a + 2i:py + 2a = (px + 2(py = ttM. (178) 

Therefore dl,dy,(fl,(py is a point on another fix-line 
identified by d%,dy. 

The initial perturbation with respect to this new point 
belonging to another fix-line reads 


A 


where now A'^^^ A'^y, A'y^, A'yy,\', lo' are evaluated at 
dl,dy, and the constants A',^,D'y are computed from 
the initial condition {Sdx)i, {Sdy)i, {S(fx)i, iS(py)i. We 
find 

{5dx)i =- ^-^sin(^'): 

LU' 


(182) 


a- ^ 2 A'aI' . 

(Aaj,)i =-^sin(^ ), 

OJ 

iS^x)i =X'{Alx + 2Aly)^cos{a - 2D'y, 

{5(py)i =)\{Aly + 2Ayy)-^ COS(^') H ^ . 


Now we use the definition of {Sdx)i, {Sdy)i, 
{S(fx)i, iS(py)i and get 

{Sdx)i = - — sin(^), 

LU 

(c~ \ 2AA , 

[day)i = -sm(^). 


{6dx)i-dx + {6dx)o al — (Sdx)o ^^(Axy + 2Ayy)Cp, =x(Axx + 2Axy)^ cos(^) - 2Dy + 2a, 


(183) 


ui 

{Sdy^i —dy (ddx)o fiy ~ (^dx)o A ^j^(Axx A 2Ayx^Cp, 

iS(fx)i =dpx A iS(px)o - v>l = {5(px)o A 2a, 

{^Ay)l =Ay A ((5(Pa;)o “ = {^Ax)o “ «• 

(179) 

We use the index 1 to denote the initial condition of the 
original perturbation now re-written with respect to the 
new fix-line. Therefore 


a;2 

(^7^y)i —^{^Axy A 2Ayy') ^ cos(^) -l- Dy a. 


Substituting Eq. 183 into Eq. 182 we find that 


C = 2{5dx)i - {5dy)i = 0 


(180) 


This means that the initial perturbation with respect to 
the point belonging to the new fix-line has the correspon¬ 
dent Cp equal to zero, and the evolution of the initial 


^ ^ A'aI' . , 

-sm(^) =-^sin(^ ), 

LU LU' 

^{Axx A 2Axy)-^ cos(^) — 2Dy 2a = 

= ^'i^'xx + ^^xy)'^ C0S(^') — 2_Dy, 
^{_-^Xy “t“ 2y4yy) — Cos(^) ~h Dy Of — 

= 2^'(A' + 2Ay )—— cos(^') -I- D' 


(184) 
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The parameters on the L.H.S. are determined by the per¬ 
turbation with respect to the original fix-line point. The 
parameters on the R.H.S are the new parame¬ 

ters that identify the original perturbation with respect 
to the new fix-line point. Multiplying the last by 2 and 
summing with the second equation we find 


1) The range of is ^ < ax < oo, because Uy > 0, 
and at da; = ^ the function ^ is not defined. 

2) For Gx +00 we always find {dx) 0. 

3) Given e a positive, very small number, we find that 
for dx = ^ + e the following approximation holds 


>lcos(^) = A'cos(^'), (185) 


which combined with the first of Eq. |184| yields 


a4'cos(^') =Acos(^) 

\t / 

A'sin{e)=^Asin{0. 

X'oj 


(186) 


From this system we obtain A', . 

Now we use the free parameter a to require that D'^ = 

0. We take the second equation of Eq. |184| and we find 

A A’ 

a = X(Axy-\-2Ayy)-^ cos{^)-\-Dy—X {A^y-\-2Ayy)-—^ cos(^ ) 

"" (187) 

substituting this value into the first equation of Eq. |184| 
we find 


-2D' = A cos(C) - A' cos(^') = 0. (188) 


187 


yields D' = 

184] 


Therefore the parameter a defined in Eq. 

0 in both second and third equations of Eq. 

We conclude that the original perturbation around the 
new fix-line point identified by a with the value in Eq. |187| 
evolves according to 


5dx 

Sdy 

5(px 

5(fy 


UJ 

9\' A' 

= -^sin(cc'0 + n, 

UJ 

cos(a;'0 -I- ^'), 

" a; ^ 

= X!{A'yx + 2 A'yy)-^ C O S ( W ^ ^ ^'). 


(189) 


That is, the initial perturbation is oscillating around the 
point d%,dy, (p% ^ o f the new fix-line identified by d^, d^. 
In fact, in Eq. |169| we simply replace all constant with 
the primed ones and set Cp = 0. This results let’s us 
proceed with the most simple case, i.e. with Cp = 0. 


XIV. APPENDIX B: PROPERTIES OF THE 
LEVEL LINES 




e = - 




(190) 


Therefore the asymptotic behaviour of the curve ^ 

" > 

< 


for dx —>■ changes according whether I ^ We 


find 



(191) 


If / = ^ then '3^ takes the special form 


'3^{dx) = 



(192) 


and for dx —>■ C"*" find Therefore 

the level line for / = ^ separates the level lines into 
two classes of curves: one class where the curves 
diverge to -l-oo when dx —>■ C"*": while in the other 
class the curves diverges to —oo when dx —>■ 

This level line is shown by the red line in Fig. 
Other two curves for / < ^ and I > ^ show the 
divergence property of 

4) If da; = / and I ^ ^ we always find 3^ = 0. There¬ 
fore for I > ^ the level curve identified by I will 
always cross 3^ = 0. 

5) The curves identified by the invariant I can have: 
1) no maximum and no minimum; 2) one maxi¬ 
mum; 3) both maximum and minimum. 

The solution of ^'(d*) = 0 for a level line defined 
by I satisfies the relation 


/ = 




3aS -C 

and the value of 3^ (d*) is given by 


6) We observe that 


^{£. + e)3r{dl) 


< 0 for / > ^ 
> 0 for i < ^ 


(193) 


(194) 


(195) 


We address here the properties of the level lines as 
defined by Eq. We distinguish the following cases. 

A. C > 0 

The level lines have the following properties: 


That means that for / > ^ the point 3^{d%) is al¬ 
ways on the opposite plane of where 3^ diverges. 
This situation is shown in Fig. The level line 
below the red curve has a maximum in (d*)+, of 
value ^((d*) + ). The point (d*)+, ^((dj)+), is lo¬ 
cated in the half plane 3^ > 0. In the other half 
the curve diverges for dx —> 
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7) The points of maximum or minimum are obtained 
by the equation = 0, which reads 

dl + {^- 3i)dx + k = 0, (196) 

and the values of the invariant for which the level 
curves can have a maximum or minimum are 

C > 9/ or ^ < 7. (197) 


8) If I is in the ranges of Eqs. 197 the two solutions 
are 


of Eq. 196 


m± = 


-^ + 37±V(C-37)2-4/C 


(198) 


We next have to control which of these solutions 
is consistent with the point 1). That defines the 
criteria if the solution is acceptable or not, (i.e. the 
criteria Eq. [^. We accept only the solution sat¬ 
isfying ^ < dx < oo. Acc ordin g to the values of 
I, and ^ the solution of Eq. 
properties 


198 


have the following 


I < 






(d*x)+ < e 

not acceptable 

(dx)- < 0 

not acceptable 

(d*x)+ > e 

acceptable 

(«:)+ = C 

not acceptable 

(d*x)- < C 

not acceptable 


(199) 


( 200 ) 


We conclude that the level lines in the region 7 > ^ 
have a maximum at (a*)+ and a minimum at ha, —> oo. 
Eor 7 < ^ the level curves have no maximum. The general 
properties of the level lines for ^ > 0 are summarized in 
the Eig. The red curve separates all the level lines 
into two classes: one class of lines without any maximum, 
and a second class of lines with a maximum. 


B. ^ = 0 

The level lines have the following properties: 

1) The range of Oa, is 0 < Oa,, because on Oa, = 0 the 
function is not defined. 

2) Eor dx —>■ +00 we consistently find ‘3^{dx) —>■ 0. 

3) Let’s study in proximity of dx = 0. We take e 
an arbitrarily small, positive number. For dx = e 
we find 


Wie) 


1 

g372- 


( 201 ) 



FIG. 27: Summary of the general behaviour of the level 
curves in the case ^ > 0. 


with i ^ 0. The curves defined by 7 are distin¬ 
guished into two classes of behaviour for the limit 
dx —f 0+ according to the value of 7: 


dx —^ 0 ~^ 


—)■ —oo for 7 > 0 
—)■ - 1-00 for 7 < 0 


( 202 ) 


If 7 = 0 the behaviour of W is different, in fact 
takes the form 


^(Oa;) = 



(203) 


and for dx —f O'*" we find W ^ oo. This particu¬ 
lar level line is shown in Fig. by the red curve. 
The curves below the red curve diverge to —oo for 
dx 0+, while the curves above the red curve will 
diverge to -l-oo for Oa, —>■ 0+. 

4) If da; = 7 then = 0. Therefore for 7 > 0 the level 
curve identified by 7 will definitely cross = 0. 


5) The equation of d*, Eq. 196 becomes now simply 
(ai(;)± = 3-^(1 i l-^l/^)/2, with solutions for —oo < 
7 < 00. 


6) Now we check if the solutions of Eq. 196 i.e. of the 
point 5) satisfy the acceptance criteria expressed in 
Eq. 83 which in 1) reads 0 < dx- According to the 
value of 7 we find the following cases 


7 < 0 


(d*x)+ = 0 

(d:)_ = 37 


not acceptable 
not acceptable 


(204) 


7 > 0 



acceptable 
not acceptable 


(205) 
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1 = 0 


(a*)+ = 0 not acceptable 
(a*)_ = 0 not acceptable 


(206) 



FIG. 28: Summary of the general behaviour of the level 
curves in the case 5 = 0- 


We conclude that in the region / < 0 the level curves 
have no maximum; in the interval / > 0 the level curves 
have a maximum at (dl.)+ = 3/ which value is i—?=. 

The description of this behaviour is shown in Fig. |28[ 


If / = 0 the behaviour of is different, in fact 
takes the form 





(209) 


and for t O’*" we find —>■ 0. This particular 
level line is shown in Fig.j^by the blue curve. The 
curves below the blues curved line diverges to —oo 
for ha; —> O'*", while the curves above the blue curve 
will diverge to +oo for ha; —>■ O’*'. 


4) If ha; = / then = 0. Therefore for / > 0 the level 
curve identified by / will cross '3^ = 0. 


5) The equation of h*, Eq. 
in the range 


196 


allows solutions for I 


or !</. 


( 210 ) 


6) Now we check if the solutions of Eq. 196 satisfy the 
acceptance criteria expressed in Eq. 83 which in 1) 
reads ha; > 0. According to the value of I, and ^ 
we find the following cases 


I<^ 


iK)+ < 0 

(K)- < C 


not acceptable 
not acceptable 


( 211 ) 


< J < 0 


(«:)+ > 0 
( 0 + = 0 
K)- > 0 


acceptable 
not acceptable 
acceptable 


( 212 ) 


0 < / : 


(a:)+ > 0 

(K)- < 0 


acceptable 
not acceptable 


(213) 


C. ^ < 0 

The level lines have the following properties: 

1) The range of ha; is 0 < ha;, because at ha; = 0 the 
function 3^ is not defined. 

2) For dx +00 we find 3^ (ha;) —>■ 0. 

3) Let’s study 3^ in proximity of dx = 0. We take e 
an arbitrarily small, positive number. For ha, = e 
we find 

3^{e) = (207) 

Therefore the curves defined by / are distinguished 
into two classes of behaviour for the limit dx —>■ O''" 
according to the value of /: 




3^ —)■ —oo for 
3^ —)■ +00 for 


/ > 0 
I < 0 


(208) 


FIG. 29: Summary of the general behaviour of the level 
curves in the case ^ < 0. 





























35 


We conclude that in the region / > 0 the level curves 
have only a maximum (a*)+; in the interval ^/9 < / < 0 
the level curves have a minimum (a*)_ and a maximum 
(a*)+; for I < ^/9 the level curves have no maximum, 
but a minimum at ax —^ oo. The description of this 


behaviour is shown in Fig. |29[ the red curve is the curve 
for 7 = ^/9. In this curve the maximum and minimum 
overlap at dx = —C/3, which is an inflection point with 

^ = l/y=3C. 
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